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EXECUTIVE SUMMARY 

The present study, which is the first of a series of investigations dedicated to 
specific issues of low pressure turbine (LPT) boundary layer aerodynamics, is aimed 
at providing detailed unsteady boundary layer flow information to understand the 
underlying physics of the inception, onset, and extent of the separation zone. A 
detailed experimental study on the behavior of the separation zone on the suction 
surface of a highly loaded LPT-blade under periodic unsteady wake flow is 
presented. Experimental investigations were performed on a large-scale, high- 
subsonic unsteady turbine cascade research facility with an integrated wake 
generator and test section unit. Blade Pak B geometry was used in the cascade. The 
wakes were generated by continuously moving cylindrical bars device. To account 
for a high flow deflection of LPT-cascades at design and off-design operating points, 
the entire wake generator and test section unit including the traversing system is 
designed to allow a precise angle adjustment of the cascade relative to the incoming 
flow. This is done by a hydraulic platform, which simultaneously lifts and rotates the 
wake generator and test section unit. The unit is then attached to the tunnel exit 
nozzle with an angular accuracy less than 0.05°, which is measured electronically. 
Utilizing a Reynolds number of 110,000, based on the blade suction surface length 
and the exit velocity, one steady and two different unsteady inlet flow conditions 
with the corresponding passing frequencies, wake velocities, and turbulence 
intensities are investigated. The reduced frequencies cover the entire operation of 
LP-turbines. The measurements were made using hot-wire anemometry. In addition 
to the unsteady boundary layer measurements, blade surface pressure measurements 
were performed at Re = 50,000, 75,000, 100,000, 110,000, and 125,000. For each 
Reynolds number, surface pressure measurements are carried out at one steady and 
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two periodic unsteady inlet flow conditions. Detailed unsteady boundary layer 
measurement identifies the onset and extension of the separation zone as well as its 
behavior under unsteady wake flow. The results presented in ensemble-averaged and 
contour plot fonns help to understand the physics of the separation phenomenon 
under periodic unsteady wake flow. 
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1. INTRODUCTION 


In recent years gas turbine engine aerodynamicists have focused their attention on 
improving the efficiency and performance of the low pressure turbine (LPT) component. 
Research works from industry, research centers, and academia have shown that reduction 
in blade number can be achieved without substantially sacrificing the efficiency of the LPT- 
blading. This reduction contributes to an increase in thrust/weight ratio, thus reducing the 
fuel consumption. Contrary to the high pressure turbine (HPT) stage group that operates in 
a relatively high Reynolds number environment, based on the operating conditions, the LPT 
experiences an adverse change in Reynolds number ranging from 50,000 to 250,000. Since 
the major portion of the boundary layer, particularly along the suction surface is laminar, the 
low Reynolds number in conjunction with the local adverse pressure gradient makes it 
susceptible to flow separation, thus increasing the complexity of the LPT- boundary layer 
aerodynamics. The periodic unsteady nature of the incoming flow associated with wakes, 
that originate from upstream blades, substantially influences the boundary layer 
development including the onset of the laminar separation, the extent of the separation zone, 
and its turbulent re-attachment. Of particular relevance in context of LPT-aerodynamics is 
the interaction of the wake flow with the suction surface separation zone. While the 
phenomenon of the unsteady boundary layer development and transition in the absence of 
separation zones has been the subject of intensive research, that has led to better 
understanding the transition phenomenon, comprehending the multiple effect of mutually 
interacting parameters on the LPT -boundary layer separation and their physics still requires 
more research. 

The objective of this present investigation under the supervision is to better 
understand the effect of the unsteady wake flow on the boundary layer behavior and the 
underlying physics of the inception, onset, and extent of the separation zone. Boundary 
layer measurements are taken on a LP turbine blade designed by Pratt &Whitney in a linear 
turbine cascade as described by Schobeiri et al.( 1995) and (1997). 
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NOMENCLATURE 


D in inside diameter of calibration nozzle (mm) 

D out outside diameter of calibration nozzle (mm) 

E anemometer voltage (V) 

m number of hot-wire samples taken per revolution of the wake 

generator 

N total number of samples of hot-wire data for time-averaging 

c blade chord length 


c P 

Pi 

P re f 

u 

V. 

V m 

V 

V 
v 
s 

Smd 


S 


o 


s 


R 


pressure coefficient, C p - , 

static pressure taps, i=1..48 
static pressure 

rotor circumferential velocity, belt translational velocity(m/s) 
free-stream velocity (m/s)- 
test section inlet velocity (m/s) 
instantaneous velocity (m/s) 

time-averaged velocity (m/s) 

turbulent fluctuation velocity (m/s) 

streamwise distance from the leading edge of the blade 

maximum separation bubble height at a streamwise distance from 

the leading edge of the blade 

starting point of the separation zone at a streamwise distance from 
the leading edge of the blade 

streamwise distance from the leading edge to the trailing edge of 
the blade 
rod spacing 
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t 

Tu 

<Tu> 

<V> 

<v> 

y 

V 

p 

a 

T 

A 

Q 


blade spacing 

time (s) Tu=yx\0to 

turbulence intensity, 

ensemble-averaged turbulence intensity 

ensemble-averaged mean velocity 

ensemble-averaged turbulent fluctuation velocity 

lateral distance from plate surface (mm) 

kinematic viscosity of air (nr/s) 

density of air (kg/m 3 ) 

cascade solidity, o= c/S B 

one wake passing period 

velocity ratio, V ax /U 

load coefficient of a constant mean diameter turbine stage, 

Mv u 2 +v u3 )/u 


non-dimensionalized unsteady parameter, 


Q = 


<P S R 


Re s local Reynolds number based on longitudinal distance, Re s =(Us)/v 

Re smd Reynolds number based on the conditions at separation 

H 12 shape factor, 

8 boundary layer thickness 

8 , displacement thickness 

8 2 momentum thickness 

83 energy thickness 

G) data uncertainty 
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2. LITERATURE REVIEW 


2.1 Flow Through a Multi-Stage Turbine 

In turbomachinery aerodynamics the motion of neighboring blade rows generates the 
periodic unsteady wakes. Several factors such as wake interaction, endwall region, and 
secondary losses affect the flow behavior and thus the performance of the turbine blades and 
turbomachine. Therefore, wake induced unsteady flow has a profound affect upon the 
overall aerodynamic behavior of blade rows. This unsteady flow influences the individual 
boundary layer of each blade and is investigated in detail by Arndt (1993) where he 
concluded that there are significant changes in the free stream turbulence levels and wake 
structures at each row location. 
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Figure 2. 1 : Schematic of the stator-rotor unsteady wake flow propagation through a 
multi-stage environment from Schobeiri et al. (1995). 
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The schematic shown in Figure 2.1 illustrates the findings by Schobeiri et al. (1995). The 
flow leaving the first stator blade row impinges upon the second blade row with a spatial 
frequency that is related to the blade spacing of the preceding stator row. Due to the rotation 
of the second row and the frame of reference change, a periodic unsteady turbulent flow 
results thus affecting the boundary layer behavior. 

As the flow leaves the first stage, two sets of wake structures are present which have 
been generated from the preceding rows. These wakes impinge upon the third row with 
different velocity defects, phase displacements, and turbulent structure. Once the flow leaves 
the third row, a new set of wakes are generated while the previous wakes have interacted 
through intensive mixing resulting in high background turbulence. The turbulence intensity 
increases for each proceeding row until some limit is reached somewhere downstream. 

2.2 Boundary Layer Transition in Turbomachinery 

Since the state of the boundary layer affects the loss production and the heat transfer 
characteristics, the turbine designer needs to know the information about the distribution of 
the laminar and turbulent flow regimes and where the transition starts and ends on the blade 
surfaces. While the transition to turbulent state has been extensively investigated, transition 
in highly disturbed environments such as turbomachinery flows are adequately documented. 
Boundary layer transition from laminar state to the turbulent state is a stochastical, three- 
dimensional, and unsteady process which occurs in a region where the flow is intennittently 
laminar and turbulent. In general, there are three significant modes of transition. The first 
is called “natural” transition, and this mode of transition starts with a weak instability in the 
laminar boundary layer and proceeds through a various stages of amplified perturbance to 
fully turbulent flow. Since this transition mode occurs only in low free-stream turbulence 
environments, it is of little importance in gas turbine flow fields. The second one is called 
“bypass” transition and is caused by a large disturbances in the external flow that directly 
produces a turbulent spot in the boundary layer. It has been the general perception that the 
presence of wakes causes the Tollmien-Schlichting waves to be bypassed (Mayle 1991). 
Extensive unsteady boundary layer investigations at Texas A&M University 
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Turbomachinery Perfonnance and Research Laboratory (Schobeiri and Radke 1994, 
Schobeiri et al. 1995) showed that despite the presence of wakes Tollmien-Schlichting 
waves are present. The third one is called “separated-flow” transition and it occurs in a 
separated laminar boundary layer. When a laminar boundary layer separates, transition may 
occur in the free shear layer near the surface and separation bubble length strongly depends 
on the transition process in the shear layer. This mode of transition occurs in gas turbine 
engines, particularly in the compressors and low-pressure turbines. 

A review of boundary layer development and the periodic unsteady flow studies are 
presented. 

2.3 Periodic Unsteady Boundary Layer Research 

One of the major effects of unsteady wake flow is on the boundary layer 
development and the transition from laminar to turbulent state. Realizing the importance of 
unsteady flow, there has been an extensive research on unsteady aerodynamics of 
turbomachines in the past two decades. Initial studies utilized simplified models to simulate 
the turbomachinery flow conditions. Precise knowledge on how the unsteady flow affects 
the distribution of the laminar boundary layer development is of great importance for the 
turbine designer and also to develop better transition correlations for use in computational 
fluid dynamic codes. 

An extensive unsteady boundary layer research program was initiated mid sixties at 
the Technical University Darmstadt, Germany. Fundamental studies by Pfeil and Pache 
(1977), Pfeil and Herbst (1979), Schroder (1983) examined the boundary layer transition on 
a flat plate with the effect of unsteady wake flow. They used squirrel cage type of wake 
creator, with a series of cylindrical rods, which is connected on two parallel rotating disks. 
Their studies show that the unsteady wakes generated by the cylindrical rods caused the 
boundary layer to become turbulent during their impingement on the plate and affect the 
length of the boundary layer transition region. They also pointed out that, wake-induced 
transition zones, which propagated downstream of the plate at velocity less than the wake 
passing velocity. 
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Orth ( 1 993) studied and quantified the effect of unsteady wake flow on the boundary 
layer transition along a flat plate. First, he observed an early onset of transition when the 
high turbulence level of the wake disturbs the boundary layer and leads to the formation of 
turbulent patches. Second, laminar becalmed regions are formed behind the turbulent 
patches so that brief periods of laminar flow are still observed beyond the location at which 
the steady flow boundary layer is fully turbulent. 

Orth ( 1 993) and Schobeiri and Radke ( 1 994) experimentally investigated the effects 
of periodic unsteady wake flow and pressure gradient on the boundary layer transition and 
development along the concave surface of a constant curvature plate. The measurements 
were carried out under zero and negative pressure gradients utilizing an unsteady flow 
research facility with a rotating cascade of rods positioned upstream of the curved plate. 
Based on the comprehensive experimental investigations, Schobeiri and Radke concluded 
that the wake induced unsteady flow significantly affects the boundary layer transition 
behavior. The onset of primary boundary layer and the location of its high turbulence 
intensity core was shifted towards the plate leading edge when subjected to periodic 
unsteady wakes. Their results were similar to those found by Liu and Rodi (1989, 1991) 
stating that the wake passing does not significant effect for the boundary layer development, 
but it cause the transition to occur earlier than in natural transition. They found that, 
transition region shifted towards to the leading edge. 

The effects of passing wakes have also been investigated by Dullenkopf et al. ( 1 99 1 ), 
Liu and Rodi (1991). Their results showed that an increase in wake passing frequency 
caused the transition region to shift towards the leading edge for the suction surface. Higher 
frequencies were noted by all the investigators to significantly enhance the heat transfer on 
both surfaces due to the increased turbulence intensity. 

Halstead et al. (1997), performed an experimental studies of boundary layer 
development on the suction surface of the low pressure turbine blade. They found a region 
of laminar and transitional flow on the suction surface without flow separation. They also 
showed that becalmed regions generated by the turbulent spots produced in the wake path, 
were effective in suppressing the flow separation. 
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Analyzing the velocity and turbulence structure of the impinging wakes and their 
interaction with the boundary layer, Chakka and Schobeiri (1999) developed an 
intennittency based unsteady boundary layer transition model by analyzing the velocity and 
the turbulence structure of the impinging wakes and their interaction with the boundary 
layer. The analysis revealed an universal pattern for the relative intennittency function for 
all frequencies and pressure gradient investigated. However, they considered the above 
investigations were not sufficient to draw a final conclusion with regard to an eventual 
universal character of the relative intennittency function. Further detailed investigations of 
the unsteady boundary layer on a high Reynolds number turbine cascade by Schobeiri et al. 
(1995) and its subsequent analysis (1997), (1998), and (2002) verified the universal 
character of the relative intennittency function. For this purpose, Schobeiri et al. (1995) 
utilized a conceptually different type wake generator. Fottner and his co-workers (2000), and 
Schulte and Hodson ( 1 996) used the same wake generating concept for investigations on the 
influence of unsteady wake flow on the LPT-boundary layer. 

Kaszeta, Simon, and Ashpis (2001) experimentally investigated the laminar- 
turbulent transition aspect within a channel with the side walls resembling the suction and 
pressure surfaces of a LPT-blade. 

Fottner et al. (2000) performed an experimental studies using mounted hot film 
measurement technique. The results shows that the transition moves when it is subjected to 
wake along the LP turbine blade, which significantly influence the loss generation in the 
suction side boundary layer. Thus, laminar separation is suppressed or reduced whenever 
the wake is present. 

Lou and Hounnouziadis (2000) investigated the mechanism of separation, transition, 
and reattachment, and the effect of main flow unsteadiness on a flat plate, contouring the 
opposite wall to create the appropriate pressure gradient on the test wall. By using a rotating 
flap which is mounted downstream of the test section, they were able to produce periodic 
oscillations of the main flow and identify the characteristic instability frequencies in the 
shear layer over a separation bubble. They show that free shear layer instability originates 
at the maximum vorticity and energy is transferred from the shear flow to the instability 
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waves where the instability waves are amplified downstream and trigger the transition. They 
also show that the momentum transfer from the shear flow is increased due to the turbulence 
fluctuations. This causes the reattachment of the flow. They investigated the Reynolds 
number effect on the transition region. They showed that for the higher Reynolds number 
results an earlier transition and reduction of the transition length, while the separation point 
does not change its location. 

Recent work presented by Cardamone et al. (2000) captures a separation zone in the 
steady state case at Re= 60,000, however, in unsteady case, the separation zone is not 
reproduced. The few boundary layer measurements are not comprehensive enough to 
provide any conclusive evidence for the interpretation of the boundary layer transition and 
separation processes and their direct impact on the profile loss, which is a critical parameter 
for blade design. 

Shyne et al. (2000) perfonned an experimental study at a Reynolds number of 
100,000 with different levels of free-stream turbulence. They showed that transition onset 
location and the length is strongly dependent on the free-stream turbulence. As the free- 
stream turbulence increase, the onset location and the length of the transition is decreased. 

Recently Schultz and Volino (2001) have shown that concave curvature has a 
significant effect on transition at high free-stream turbulence. 

Volino and Hultgren (2001) performed an experimental study and measured the 
detailed velocity along a flat plate which is subjected to the same dimensionless pressure 
gradient as the suction side of the modem LP turbine blade. They showed that the location 
of the boundary layer separation does not strongly depend on the Reynolds number or free- 
stream turbulence level, as long as the boundary layer remains non-turbulent before 
separation occurs. But they also showed the extent of the transition are strongly dependent 
on the Reynolds number and turbulence intensity. With higher free-stream turbulence levels, 
transition appears to occur in a bypass mode, similar to the attached boundary layers. On the 
other hand, under low turbulence intensity, transition is induced by instability waves in the 
shear layer of the separation bubble. 
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Treuren et al.(2002) performed an experimental study at a very low Reynolds 
number of 25,000 and with different free-stream turbulence levels. They showed that a 
massive separation at very low Reynolds number of 25,000 is persistent, in spite of elevated 
free stream turbulence intensity. However, at a higher Reynolds number of 50,000, there is 
a strong separation on the downstream of the suction side for the low free-stream turbulence 
level. The separation zone eliminated for the higher free-stream turbulence level of 8-9%. 
And they also added vortex generators to provide mixing for the reattachment. At very low 
Reynolds number of 25,000, they showed a massive separation on the suction surface, in 
spite of elevated free-stream turbulence level and vortex generators. However, elevated free- 
stream turbulence and vortex generators helps the flow to reattach to the surface at higher 
Reynolds number. 

To prevent the separation in LP turbine vortex generator jets and dimples can be 
used. Recently Lake et al. (2000) suggested to use dimples on the suction surface of a 
turbine blade to control separation for the low Reynolds number. Their research showed a 
decrease in loss coefficient at low Reynolds number with an aerodynamic penalty at higher 
Reynold number. Another method, suggested to control or prevent the separation is the use 
of vortex generator jets. This method produce no important adverse effects at higher 
Reynolds number. Recently Bons et al. (1999) applied this method to low pressure turbine. 
He showed that the separation region was almost eliminated and there is a significant 
reduction in momentum loss downstream of the blade. Both grooves and vortex generators 
has been applied to wind turbines. While energizing the boundary layer and preventing the 
separation at low Reynolds numbers, there is a penalty at high Reynolds number which is 
considered unacceptable. On the other hand, Lin (1999) showed that micro-vortex 
generators are effective reducing the separation zone with a little drag penalty. 

Even though, the published studies helped in better understanding the wake 
development, unsteady transition and the boundary layer development including the onset 
and extent of the laminar separation and its turbulent reattachment, they are neither complete 
to be applicable to turbomachinery flow nor to develop a theoretical framework. 
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3. OBJECTIVES 


The objective of this study is to understand the physical mechanisms associated with 
the effects of wakes on separated LP turbine flows along the blade surface and to create high 
fidelity data base suitable to be used for validation of transition, turbulence, separation 
models as well as CFD-codes. Measurements were performed in a single row cascade 
representative of LP turbine designs. A large-scale, high subsonic research facility has been 
used for the investigation. The wake passing frequency is varied from 0 to 120 Hz to 
simulate various turbomachinery-operating conditions. 

Careful documentation was made of (1) the development of the boundary layer 
which approaches the separated flow region, (2) the free shear layer over the separation 
region and (3) the reattached flow downstream of the separation bubble. 

Detailed steady and unsteady boundary layer velocity, turbulence intensity, and static 
pressure measurements have been performed along the suction and pressure surface of one 
blade. A single hot wire probe has been used for the velocity measurements. The results are 
presented in temporal-spatial domain for a detailed display of transition and further 
development of boundary layer. 

In addition to unsteady boundary layer measurements, blade surface measurements 
were investigated at Re=50, 000, 75,000, 100,000, 110, 000 and 125,000. For each Reynolds 
number, surface pressure measurements are carried out at one steady and two different 
unsteady inlet conditions. 
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4. EXPERIMENTAL TEST FACILITY 


For the detailed investigations on the varying wake passing frequencies, a large 
scale, high-subsonic research facility has been used. The facility in its original configuration 
is described by Schobeiri et al. (1995) and (1997). The facility is shown in Figure 4.1. The 
facility consists of a large centrifugal fan, straight pipe, diffuser, settling chamber, nozzle, 
cascade test section, motor driving the wake generator. A detailed description of the test 
facility is given below. 

4.1 Air Supplier 

Air is supplied from a large centrifugal fan, which is driven by 93 kW 3-phase AC 
electric motor running at a rated speed of 1760 rpm, with a volumetric flow rate of 15 m\ 
The blower reduces the rotary speed to a fan speed of 1035 rpm by a belt-pulley 
transmission. The capacity of the fan is 100 m/s for the maximum mean velocity at the test 
section. It corresponds to a Mach number of M=0.6 for a turbine cascade with an 
acceleration ratio of V 2 /V 1 =2. Various inlet velocity can be obtained by changing the fan 
speed or the throttle which is located at the exit of the fan. For the present investigation, the 
throttle is adjusted to give a maximum Reynolds number 120,000 at the inlet of the test 
section. To protect the hot-wire system, a fiber-glass filter with a thickness of 50 mm and 
a pore size of 5pm is attached at the inlet of the fan. 

Since the fan exit has a higher elevation than the test facility, a set of transition ducts 
followed by a straight pipe with a diameter of 600 mm are used for the smooth channeling 
of the flow from the fan exit to the diffuser inlet. A Prandtl probe is located at the inlet of 
the diffuser to get the variations in the mean velocity resulting from the voltage difference, 
fluctuations of the fan speed, or change of environmental temperature and pressure. 

4.2 Diffuser, Settling Chamber, and Nozzle Assembly 

The straight pipe leads to the diffuser with a length of 1500 mm and area ratio of 
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Figure 4. 1 : Turbine cascade research facility with the components and the adjustable test 
section. 

1 :4.75. The flow is decelerated before entering the settling chamber, which consists of four 
sections. Each section has a length of a 900 mm and cross section of 1500x900 mm. Five 
screens and one honeycomb flow straightener are used in the settling chamber to control the 
flow uniformity and turbulence intensity level. The screens are made of stainless steel and 
have different mesh sizes and wire diameters. The first screen is placed at the diffuser exit 
to reduce the scale of vortices generated by the diffuser vanes. The honeycomb flow 
straightener with a cell size of 6 mm is placed between the first two settling chamber. The 
second screen is located directly downstream of the honeycomb. The honeycomb and screen 
combination results in a much lower turbulence intensity than a honeycomb alone because 
the large eddies exiting from the honeycomb cells are broken into smaller eddies. The other 
screens are placed downstream of the second, third and fourth sections of the settling 


NASA/CR— 2009-2 14831 


13 




chamber. The nozzle is located downstream of the settling chamber with an area ratio of 
6.75 : 1 . The flow accelerates before entering the wake generator and the nozzle provides a 
smooth transition of flow from the settling chamber to the wake generator. A straight 
channel transition piece is attached to the nozzle to accommodate the test section 
inclination. 

4.3 Periodic Unsteady Inlet Wake Flow Generator 

The translational motion of the wake generator generates the periodic unsteady 
wakes. It consists of a series of cylindrical rods that are placed on both side of the timing 
belt that is driven by an electric motor. The timing belts with a length of 5 m span over 1 0 
pulleys arranged around the cascade test section. The drive pulleys also control the belt 
tension. For the present investigation, a certain number of cylindrical rods with a diameter 
of 2 mm are mounted between two timing belts. The purpose of the use of cylindrical rod 
is to get the similar rotor blade wakes. Attaching the rods to the belts or detaching the rods 
from the belt can change, the rod spacing, S R . The wake flow generator system is driven by 
an electric motor with a maximum power of 7.5 kW (lOhp). The controller allows the belt 
translational speed up to 25 m/s. For the present investigation, a belt speed of U=5 m/s is 
chosen. And also, a fiber optic system is used to monitor the wake passing frequency and 
it also serves as a trigger for the data acquisition. 

4.4 Sequential Generation of Unsteady Wakes 

With the special design of the test facility, it is possible to measure the velocities for 
different wake passing frequencies in one data set, when perfonning the boundary layer 
experiments. The wake passing frequency can be varied by changing the speed of the belt 
or the spacing of the rods on the belt, S R , shown in the Figure 4.2. For the present study, the 
belt is subdivided into three different clusters of the same length. The first one with the 
spacing S R =°°, which is steady flow case with no rods. The second and the third part have 
the rod spacing of S R = 160mm and 80 mm respectively. Within each cluster the rods have 
the same spacings. The clusters are all arranged with a certain distance from each other. By 


NASA/CR— 2009-2 14831 


14 



using the triggering system, a continuous data set is acquired and the buffer zones between 
the data clusters are visible. The data analysis program is able to cut the buffer zone and 
evaluates the data for each cluster. 



Fig. 4.2: (left)Wake generator with different rod spacing clusters, (right) velocities 
produced for rod spacing of S R = 80 mm (top), 160 mm (middle), and °°. 


In order to account for the unsteadiness caused by the frequency of each individual 
wake generating cluster and its spacing, the flow velocity and the cascade parameters the 
unsteady flow parameter Q is defined: 
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Many researchers have used Strouhal number as an unsteady flow parameter which 
includes the wake generator speed and the inlet flow velocity. However the unsteady flow 
parameter which is defined above and includes the cascade solidity, the flow coefficient, the 
blade spacing and the rod spacing, is an extension of Strouhal number in the sense that it 
additionally includes the rod spacing S B and blade spacing S R . The individual cluster 
configurations with the corresponding Q parameter are specified in Table 4.1. 

Table 4.1: Test section cascade parameters 


Parameters 

Values 

Parameters 

Values 

Inlet velocity 

V in = 4 m/s 

Inlet turbulence intensity 

Tu in = 1.0% 

Rod translational speed 

U = 5.0 m/s 

Blade Re-number 

Re = 110,000 

Nozzle width 

W = 200.0 mm 

Blade height 

h B = 200 mm 

Blade chord 

c = 203.44 mm 

Blade stagger angle 

y = 26° 

Blade axial chord 

c M = 182.85mm 

Cascade solidity 

Q= 1.248 

Blade suction surface length 

L ss = 270.32 mm 

Cascade angle 

(p = 55° 

Cascade flow coefficient 

$ = 0.80 

Cascade spacing 

S B = 163 mm 

Inlet air angle, Fig. 2 

a, = 55.0° 

Exit air angle, Fig. 2 

a 2 = 145.0° 

Rod diameter 

D r = 2.0 mm 

Rod distance to lead, edge 

L r = 122 mm 

Cluster 1 (no rod, steady) 

S R = 00 mm 

Q - parameter steady case 

II 

o 

o 

Cluster 2 rod spacing 

S R = 160.0 mm 

Q - parameter for cluster 1 

Q= 1.59 




•f?— 3.18 


4.5 Cascade Test Section 

The cascade test section shown in Figure 4.3 is located downstream of the wake 
generator and it has an inlet cross section of 1000 mmx200 mm. The test section can 
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incorporate to 5 blades with a blade height of 200 mm and a chord of 203.44 mm. For 
boundary layer investigations, five identical blades designed by Pratt& Whitley were 
implemented whose cascade geometry is given in Table 4.1. The blade resembles the 
essential feature such as the laminar boundary layer separation that is inherent to typical 
LPT-blades. This blade number is necessary and sufficient to secure a spatial periodicity for 
the cascade flow. It also consists of two plexiglas side walls. One of the side walls includes 
inlet, exit traversing slots and the slots for the boundary layer measurements which are 
placed between the second and the third blade from the bottom as shown in Figure 4. 1 item 
16. 




Figure 4.3: Cascade geometry, flow, and stagger angle are listed in Table 1. Number 
of blades=5, instrumented blades are: SPB-1 and SPB-2 with static pressure taps, 
HFB: Blade with surface mounted hot films to be used for future investigations 


4.6 Probe Traversing System 

A computer controlled traversing system is used for all the inlet, exit and boundary 
layer measurement at the cascade test section. The traversing system is mounted vertically 
on one of the plexiglas side walls. It consists of a power screw with a maximum traversing 
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length of 1200 mm which is connected to a stepper motor to provide the drive. It receives 
a signal from the LAB VIEW program and turns the power screw by an exact number of 
pulses required. An optical encoder with 400 radial slots connected to the stepper motor 
provides a feedback to the stepper motor for accurate movement along the traversing system 
and the traversing system is capable to move in a small steps of 2.5 |lm. The traversing 
system also includes a switch to prevent the probes movements beyond the desired limit. 

4.7 Hydraulic Platform 

To account for a high flow deflection of LPT-cascade, the entire wake generator and 
test section unit including the traversing system was modified to allow a precise angle 
adjustment of the cascade relative to the incoming flow. This is done by a hydraulic 
platform, which simultaneously lifts and rotates the wake generator and test section unit. 
The unit is then attached to the tunnel exit nozzle with an angular accuracy better than 0.05°, 
which is measured electronically. 


NASA/CR— 2009-2 14831 


18 



5. INSTRUMENTATION AND DATA ACQUISITION 


The test facility is instrumented for fully automated data acquisition and is similar 
to the data acquisition system described in Schobeiri, et al. (1995). 

5.1 Instrumentation of the Test Facility 

The instrumentation and data acquisition concept has been further developed and 
applied to this research facility. The schematic layout of the test facility instrumentation is 
shown in Figure 5 . 1 . For the aerodynamic investigations, the data acquisition and reduction 
system is controlled by a personal computer with a Pentium 3 processor. For the data 
acquisition and analysis, 16 channel 12-bit analog-digital (A/D) board (NI PCI-MI O 16-El) 
is installed in the computer. The outputs from instruments are connected to the channels of 
the NI extension board from where the data are sampled by the computer. The NI board also 
has a separate DAC channel through which an analog voltage can be sent out. This DAC 
channel has been used in a sequential data acquisition from a 48-port scanivalve system for 
the pressure measurement. 

Mean velocities and turbulent fluctuations are obtained using two existing 
commercial 3-channel (TSI, IF A- 1 00), constant temperature hot-wire anemometer system. 
The systems have signal conditioner with an adjustable low pass filter, DC-offset, and 
adjustable gain. The offset and the gain can be adjusted to get the maximum resolution and 
signal noise ratio. The low pass filter of the signal conditioner has been set to the half of the 
sampling frequency of 20 kHz to avoid the aliasing effects. The hot wire has 250° C sensor 
operating temperature. All the hot wire sensors used are made of tungsten with a diameter 
of 4 pun. The first channel of the first anemometer is connected to a single hot-film probe 
for measuring the time-dependent velocity distributions within the blade channel. The 
second channel is used for calibration of the single hot-wire probe connected the second 
channel of the Ni extension board. 

To measure the free-stream temperature, an OMEGA T-type thennocouple mounted 
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at the top of the Plexiglas wall is connected to the third channel of the Nl-board. Fourth 
channel of the Nl-board is connected to a high precision capacitance type 0-10 ininHg 
differential pressure transducer (MKS) with a 0-10 Volt analog output. A Prandtl probe, 
placed upstream of the diffuser, monitors the reference velocity at a fixed location. 

For the unsteady investigations, a trigger is required for phase-averaging the 
instantaneous data. A fiber optic proximity sensor placed close to the pulley of the wake is 
used. The passage signals of the rods are detected by the sensor using a silver-coated 
reflective paint on one side of the wake generating timing belts. This sensor gives an 
accurate readout of the speed of the wake generator and the passing frequency of the rods. 

For static pressure measurement and also to make measurements with the five hole 
probe, it is necessarily to read from several ports. A 48-port scanivalve system has been used 
for this purpose. The output of the scanivalve system is connected to the MKS pressure 
transducer. By using the scanivalve system, it is possible to measure 48 pressures using the 
same pressure transducer. 

5.2 Data Acquisition Procedure 

The data acquisition system is computer controlled. The required code is written 
with LAB VIEW. The computer acquires the data as a voltage from the differential pressure 
transducer, thermocouple, trigger circuit and from TSI IFA-100 flow analyzer. 

For unsteady boundary layer investigations, 20,000 samples were taken at a rate of 
20 kHz for each of the 100 revolutions of the wake generator. The low pass filter of the hot 
wire anemometer is set to 10 kHz to prevent the aliasing effect. 

One of the blades was produced with 48 static pressure taps to obtain pressure 
distribution on the suction and the pressure surface of the blade. The 48 taps on the blade 
are connected with the 48 ports of the scanivalve and the output is connected to a high 
precision differential pressure transducer. When the data acquisition program is executed, 
the computer acquires 4096 samples at a rate of 4096 Hz from the first port and calculates 
the mean pressure. Then the program sends out a voltage pulse through DAC port and steps 
the scanivalve to the next port. The system waits until steady state is reached before 
sampling the data of the next port. The process is repeated until all the ports are scanned. 
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6. CALIBRATION PROCEDURE OF THE SINGLE HOT WIRE 


Calibration of sensors is of great importance for accurate measurement of any 
physical quantity. The data acquisition procedures involve sampling of voltage signals from 
the NI board channels. These voltage signals have to be converted into real physical 
quantities by means of calibration. These quantities are temperature, pressure, flow velocity, 
flow angle etc. The pressure transducer is calibrated using a dead weight tester while the 
thennocouple is connected to the thermocouple calibrator. Least square polynomial fits are 
used to fit the averaged output voltages with the actual readings obtained from the dead 
weight tester and thermocouple calibrator. These coefficients are stored in separate data 
files. However, the calibration of five-hole probes and hot wire probes are more extensive 
and is perfonned in unifonn flow open jet calibrating facility. 

6.1 Description of the Calibration Facility 

A schematic of the calibration facility is shown in Figure 6. 1 . This facility has been 
used to calibrate the hot wire and five-hole probes. Compressed air is drawn from a reservoir 
passes through a pressure regulator, filter and flow control valve before entering the 
calibration facility. It consists of settling chamber followed by a pipe with three axis- 
symmetric sections, each having a diameter of 1 50 mm. Screens are placed between flanges 
of the sections to reduce the turbulence level. The turbulence intensity level at the jet exit 
is approximately 0.35%. A nozzle with an exit diameter of 38. 1 mm is attached at the end 
of the pipe. The inlet and outlet of the nozzle are parallel to the axis of the facility. A 
thennocouple is placed in the first section of the calibration facility to measure the air 
temperature. And a pressure tap is placed in the last section of the pipe to measure the 
pressure difference between static pressure inside the pipe and the atmospheric pressure with 
a differential pressure transducer. Two lynx stepper motor controllers are used to 
automatically change the pitch and yaw angle. The tip of the probes were placed about 0.4 
diameters downstream of the nozzle exit plane and within the potential core of the jet. The 
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Figure 6.1: Test Facility for calibrating hot-wire and the five-hole probes 


analog signal from the differential pressure transducer, thermocouple, stepper motor control 
and hot wire anemometer are transferred to the A/D board from where they are sampled. 

6.2 Description of Calibration Technique for Single-Wire Probes 
6.2.1. Single Sensor Hot -Wire Probe 

To measure the unsteady velocity in the boundary layer, a single sensor hot wire 
probe is used which is shown in Figure 6.2. The wire is made of tungsten with a diameter 
of 4 pm. The wire temperature is constant, and the sensor temperature of the T probe =250° C. 
The amount of electrical energy dissipated in the sensors is a measure of the cooling effect 
of the fluid past the heated sensor. The actual velocity is a function of the dissipated energy 
and the fluid temperature, assuming constant pressure and fluid properties. 

A constant temperature anemometer is used to vary the electrical current through the 
hot wire in order to maintain a constant wire temperature T probe . This is done by controlling 
the hot wire’s resistance, which is a linear function of the temperature. 
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Figure 6.2: Single sensor hot-wire probe 


6.2.2 Calibration Procedure 

During the velocity calibration, single hot wire probe is kept nonnal to flow and the 
flow velocity is varied. The exit velocity is calculated by the differences between the static 
pressure at the inlet pj and the static pressure at the outlet of the nozzle p 2 . From the 
Bernoulli-equation: 

Pi ~ Pi = y0 2 2 - Vi 2 ) (6.1) 


and the continuity equation for incompressible flow: 


v x d x 2 - Vr,d 2 


the formula for the exit velocity is: 


\2Ap d 4 
p d 4 -d 2 4 


the density is calculated by the ideal gas equation: 


The temperature is measured by a thermocouple inside of the facility. All variables 
are displayed in Figure 6.3. 

The exit velocity of the nozzle was varied from 0 to 12 m/s. For each velocity, a set 
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of data from constant temperature anemometry, thermocouple, pressure voltage signals 
were taken. The data are averaged for each channel. The correlation between the hot wire 
voltage (E) and the velocity (V) by a fourth-order polynomial: 

V — Uq + E + CI2E + u^E + a 4 E (6.5) 

The coefficients are calculated using the least-square method. The results of the calibration 
is shown in Figure 6.4. 



Figure 6.3: Nozzle of calibration test facility 



Figure 6.4: Calibration curve of single hot-wire probe 
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The temperature compensation of the velocity is made later, adjusting the anemometer 
output voltage according to the formula. 



= E 


T - T 

probe cal 

T - T 

probe 


(6.6) 


where E is the sampled anemometer voltage, T is the actual temperature and T probe is the 
operating temperature of the hot wire. 
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7. DATA REDUCTION 


7.1 Data Reduction for Boundary Layer Measurement 

7.1.1 Steady Flow 

For data reduction and analysis, the characteristic response of the hot wire probe is 
stored in form of calibration constants. The instantaneous velocity component is calculated 
from the temperature compensated instantaneous voltages by using the calibration 
coefficients. The instantaneous velocity is represented in the following form: 

V=V + v (7-1) 

where V is the mean (time-averaged) velocity and v is the turbulent fluctuation component. 
The mean velocity, also known as the time-average, is given by: 

V ={ jt V J P.2) 

V j=i 

where N is the total number of samples taken at one boundary layer location. A sampling 
rate of 20kHz was used for investigating the steady flow (no wakes). The root mean square 
value of the turbulent velocity fluctuation is obtained from the instantaneous and mean 
velocities by: 

v= M (r '~ Fy <73) 

and the local turbulence intensity is defined as: 

v i rr^ 37 

Tu, = -=x 100 = ■= —V (V, - V) x 100 (7.4) 

V V V N J 
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7.1.2 Periodic Unsteady Flow 

The unsteady data were reduced by ensemble-averaging method with respect to the 
rotation of the wake generator. At each boundary location, 20,000 samples were taken at a 
sampling rate of 20 kHz. The ensemble-averaged results were calculated for 100 rotation 
of the wake generator. For unsteady cases, the ensemble-averaged velocity, fluctuation 
velocity, and the turbulence intensity were calculated from the instantaneous velocity 
samples by: 

( 7 . 5 ) 

7=1 

v, (',-)*< v, (/,)>= rjrlXU, )- < v,(t,)>] 2 (7.6) 

V N M 


Tu, (t, ) =< Tu j (t i ) >= 


<V i {t i )> 


x 100 


(7.7) 


where N= 100 is the total number of wake generator cycles and M the number of samples 
taken per cycle. <V I (t,)> is the ensemble averaged velocity for a particular spatial location. 
The ensemble-averaged boundary layer parameters such as displacement thickness <Sj>, 
momentum thickness <S 2 >, and shape factor <H 12 > are calculated from: 


<s> i 


S 1 =< >= f 1- 


< v> 

< V >, 


dy 


(7.8) 


< 5 > 


S 2 =< S 2 >= J 


<v> 


<v> 


<v> 


<v> 


dy 


(7.9) 


H\2 >- 


< 8 X > 
c 8 2 > 


(7.10) 
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7.1.3 Accurate Estimation of Boundary Layer Thickness 

The accuracy of the boundary layer parameters such as displacement, momentum, 
and energy thickness depends primarily on an accurate estimation of the boundary layer 
thickness S. Since in a turbine or compressor blade channel the velocity distribution outside 
the boundary layer is not constant as shown in Figure 7.1 (a). The estimation method used 
in conventional flat plate boundary layer calculation leads to incorrect results. In the 
following, a simple iterative method, developed by Schobeiri (2003), is presented that 
accurately detennines the boundary layer thickness and thus the boundary layer parameters. 



Figure 7.1: (a): Boundary layer velocity traverse, (b) Boundary layer edge velocity on 
pressure and suction side from Schobeiri (2003) 


In Figure 7.1 the experimental velocity distribution on the pressure and suction 
surfaces are shown. To determine the boundary layer thickness, first a least square fit (LSQ) 
is performed and the velocity profile outside the boundary layer and its intersection with the 
u-axis is found. Starting from a first guessed value 5 1 , the area F 1 is found using the relation : 

s' 

f' = \ (u,„, - u )dy urn 

0 

In the next iteration step with §' the boundary of the integral in (7.1 1) is increased 
by a very small increment (for example 0.001mm) and the area F is calculated: 


NASA/CR— 2009-2 14831 


29 



(7.12) 


S' 

F ‘ = I (W< - «)<fy 

o 

We define the following ratio with the accuracy C. 

\F' -F m \ 

R= 7 L < e= 10' 6 

F 


(7.13) 


and check if the € has been reached. If the next iteration step with 6 1+1 fulfils the above 
criterion, then the desired boundary layer thickness is § n , which can be used or the boundary 
layer thickness to be implemented into equations (7.8) to (7.10). 
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8.UNCERTAINTY ANALYSIS 


Uncertainty analysis is performed on the single hot-wire measurements after 
calibration and data reduction utilizing the methods of Kline-McClintock (1956) and S. 
Yavuzkurt (1984). Kline and McClintock presented a technique to calculate how the 
uncertainty in individually measured values add together to result in an uncertainty for a 
quantity calculated from those values. However, it should be noted that the determined 
uncertainty may not be the actual uncertainty of the results because the Kline and 
McClintock method calculates the worst possible uncertainty. The main difference between 
these two uncertainty analysis method is coming from the way in which the actual 
uncertainty is calculated. The uncertainty for the Kline and McClintock method is a 
equivalent to the 95% confidence band that surrounds the single hot-wire curve fit. Since 
each point may have a different confidence interval, the uncertainty for each data point may 
be different. On the other hand, S. Y avuzkurt calculates the uncertainty by taking the square 
root sum of the squared errors divided by the number of data points. The results of the 
uncertainty analysis and the most significant equations used to determine the uncertainty for 
each case is presented below with the actual calculations documented in Appendix B. 

8.1 Uncertainty Analysis for Single Hot-Wire Measurements 

The uncertainty in velocity for the single hot-wire probe after calibration and data 
reduction is given in Table 8.1 for the Kline and McClintock and Table 8.2 for the S. 
Yavuzkurt method. The uncertainty in velocity during calibration increases when the flow 
velocity decreases. It is primarily due to large uncertainty in the pressure transducer that is 
used for calibration. As for the uncertainty in velocity after data reduction, the trend 
continues due to the fact that the calibration uncertainty is incorporated into the final 
velocity. For this case, the pneumatic velocity uncertainty and the total calibration 
uncertainty was determined using equations 8.1 and 8.3. 
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The uncertainty in the pneumatic velocity was found with the following equation; 


CO., = 


v dv^ 

K dP J 


CO , 


+ 


f dv^ 

v dp J 


1 2 


CO. 


+ 


< dv} 

dD„ 


CO 


D„I 


+ 


'avp 

,)D. 


CO 


D,„ 


( 8 . 1 ) 


In the above equation, G) p is the uncertainty of pressure readings from pressure transducer, 
tO p is the uncertainty during the calculation of density, 0) Dno 0) Dni are respectively the 
uncertainty of the outer and inner diameter of calibration nozzle. The pneumatic velocity is 
calculated: 


v = 


IP 




f D ^ 

4 

f 

1- 

v D ni ) 



( 8 . 2 ) 


In the above equation, P is the transducer pressure, D no is the outside diameter of the 
calibration nozzle, D ni is the inside diameter of the nozzle, p is the air density. The 
uncertainty in the effective velocity after the calibration for the Kline and McClintock 
method is: 


®Ve, ~ 


dVe 

oe 1 


\ 2 


+ CO cf 2 + COy 1 


(8.3) 


where (0 Ve is the effective velocity, G) e is the uncertainty of the constant hot wire 
anemometry voltage readings, (O v is the velocity and (O cl is the uncertainty of the least 
square fit and equivalent to the 95% confidence interval which surrounds the curve fit. The 
effective velocity for any particular wire is represented by the following equation: 

Vc’j = a 0j + a Xj E . + a 2j E 2 + a 3j E 2 + t/ (8.4) 


where a ;j are the coefficients determined by the least square fit. 

As for the uncertainty in velocity after data reduction, it is given by the following equation; 
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1 de 


\ 2 


CO 


+ CO 


calib 


where G) calib is the uncertainty in the calibration. 

The uncertainty in the curve fit from S. Yavuzkurt method is given by; 


c,= 


I 

i = 1 


AV 

Ve) 


\ 2 


*i 


n 


(8.5) 


(8.6) 


where n is the number of data points, and AV is the residual. This value is squared and 
added to the square of AV/V, and raised to the one half power. This calculated value is the 
uncertainty of effective velocity. 

The curve fit uncertainty for the S. Yavuzkurt method is much larger than the 
uncertainty determined from the 95% confidence interval. 

Table 8.1: Uncertainty in velocity after calibration and data reduction 


(Kline and McClintock Method) 


V (m/s) 

GJV/V (%) 

«V e /V e (%) callb 

WV/V (%) red 

3 

5.59 

5.92 

6.07 

5 

2.06 

2.36 

2.44 

10 

0.54 

0.76 

0.83 


Table 8.2: Uncertainty in velocity after calibration and data reduction 

(S. Yavuzkurt Method) 


V (m/s) 

G)V/V (%) 

«V f /Ve(%)cal lb 

WV/V (%) red 

3 

5.59 

5.71 

5.72 

5 

2.06 

2.30 

2.35 

10 

0.54 

1.28 

1.37 
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9.RESULTS AND DISCUSSIONS 


In this section, the experimental results from the steady and unsteady boundary layer 
development, periodic unsteady wake development, surface pressure measurements and the 
velocity distribution will be presented at different Re-numbers. 

Detailed steady and steady boundary layer and surface pressure measurements were 
taken at four different Reynolds number to investigate the influence of periodic unsteady 
wake flow on the boundary layer development along the suction and the pressure surface of 
the LP turbine blade and the onset of the separation zone. For each Reynolds number three 
different frequencies were performed. For the generation of the unsteady wake flow along 
the pressure and the suction surfaces of the LP turbine blade, two clusters of cylindrical rods 
with the diameter d R =2mm were attached to the two belts of the wake generator. The 
cylindrical rod diameter is chosen according to fulfill the generation of the drag coefficient 
C D which is approximately equal to the drag coefficient of the turbine blade with the chord 
and spacing given in the Table 4.1. 

To account for the unsteadiness caused by the frequency of the individual wakes, the 
flow velocity, the spacings and the cascade parameters, an unsteady flow parameter or 
reduced frequency Q is defined which is including the cascade solidity a, the flow 
coefficient (j), the blade spacing S B , and the rod spacings S R . 


c U t7 ^ 


c 

with <T = and (p 

S B 


V 

ax 

U 


Experimental investigations were performed for three different reduced frequencies 
Q=0.0, 1.59, and 3.18 that correspond to the rod spacings S R =°° mm, 80 mm, and 160 mm. 
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9.1 Surface Pressure Distributions 

Detailed pressure measurements were taken by a 48-port scanivalve at four different 
Re-numbers of 50,000, 75,000, 100,000, 1 10,000 and 125,000. For each Reynolds number 
three different reduced frequencies of Q=0, 1.59, 3.18 were applied that correspond to rod 
spacings of S R =oo ,160 mm, and 80 mm. The pressure signals inherently signify the time- 
averaged pressure because of the internal pneumatic damping effect of the connecting pipes. 
The differences in pressure distributions between steady and unsteady cases, on the suction 
surface, are caused by the drag forces of the wake generated by moving rods. The drag 
forces cause a momentum deficiency that lead to a reduction of total and static pressure. 

The time-averaged pressure coefficients along the pressure and suction surfaces are 
plotted in Figure 9. 1 and 9.2. Starting from Re=50,000 and steady state with Q=0, shown 
in Figure 9.1 (a), the suction surface (upper portion), exhibits a strong negative pressure 
gradient. The flow accelerates at a relatively steep rate and reaches its maximum surface 
velocity that corresponds to the minimum C p =-4.0 at s/s o =0.42. Passing through the 
minimum pressure, the fluid particle within the boundary layer encounter a positive pressure 
gradient that causes a sharp deceleration until s/s o =0.55 has been reached. This point 
signifies the begin of the laminar boundary layer separation and the onset of a separation 
zone. As seen in the subsequent boundary layer discussion, the separation zone 
characterized by a constant c p - plateau extends up to s/s o =0.746, thus occupying more than 
21% of the suction surface and constituting a massive separation. Passing the plateau, the 
flow first experiences a second sharp deceleration indicative of a process of re-attachment 
followed by a further deceleration at a moderate rate. On the pressure surface, the flow 
accelerates at a very slow rate, reaches a minimum pressure coefficient at s/s o =0.42 and 
continues to accelerate until the trailing edge has been reached. Unlike the suction surface, 
the pressure surface boundary layer does not encounter any adverse positive pressure 
gradient, therefore it remains fully laminar. 

Considering the unsteady case with the flow reduced frequency Q=1.59 
corresponding to a rod spacing of S R =160 mm. Figure 9.1 (a) shows two important and 
different characteristics : ( 1 ) There is a significant difference in pressure distribution between 
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steady and unsteady cases. As mentioned above, this deviation is attributed to the 



0 0.2 0.4 0.6 0.8 1 

s/s 0 



0.2 0.4 0.6 0.8 

s/s 0 


Figure 9.1: Static pressure distribution at (c) Re=100,000 and (d) Re=l 10,000 and 
reduced frequencies Q=0 (S R =°°), Q=1.59 (S R =160 mm), Q=3.18 (S R =80 mm) 

momentum deficiency that leads to a reduction of total and static pressure. (2) For 
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Figure 9.2: Static pressure distribution at (e) Re=125,000 and reduced 
frequencies Q=0 (S R =°°), Q=1.59 (S R =160 mm), Q=3.18 (S R =80 mm) 


Re=50,000 and less than 125,000, the wakes have a substantial reducing impact of the 
separation zone. Increasing the reduced frequency to Q=3. 18 by reducing the rod spacing 
to S r =80 mm causes a slight shift of the C p -distribution compared with other frequencies. 

Increasing the Reynolds number to Re=75,000, has not brought major changes in 
steady state C p -distribution. This is also true for the subsequent higher Reynolds number 
cases at steady state Figure 9. 1 and 9.2 (b,c,d,e). However, the higher Reynolds number with 
unsteady wakes reveals that the noticeable deviation in pressure distribution between the 
steady and unsteady cases discussed above is diminishing by increasing the Re-number as 
shown in Figure 9.1 and 9.2 (b,c,d,e). There is two factors which is contributing to this 
deviation: (1) As mentioned above, the momentum deficiency and the associated pressure 
losses caused by moving wakes. (2) The boundary layer energizing effect of the impinging 
wakes. Although the impinging wakes cause velocity and momentum deficits, their high 
turbulence intensity vortical cores provide an intensive exchange and transfer of mass, 
momentum, and energy to the blade surface, thus energizing the low energetic boundary 
layer. In conjunction with the surface pressure distribution, the kinetic energy of the lateral 
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velocity fluctuation plays a crucial role. In case of a low Re-number flow, the strong 
damping effect of the wall shear stress has the tendency to reduce the lateral turbulence 
kinetic energy, thereby diminishing its surface pressure augmenting effect. Increasing the 
Re-number results in a decrease of the damping effect of the wall shear stress, allowing the 
kinetic energy of the lateral velocity fluctuation component to be added integrally to the 
surface pressure, thus offsetting the wake deficit effects on the pressure distribution. This 
fact is clearly shown in Figure 9.1 and 9.2 (a,b,c,d,e), where the pressure distributions of 
unsteady flow cases at Q=1.59 and Q=3.18, with the exception of the separation zone, 
systematically approach the steady state cases at Re=75,000, 100,000 and very visibly at 
Re=125,000. It is worth noting that the impact of the unsteady wakes on the extent of the 
separation zone is preserved regardless of the Reynolds number variation performed in this 
study. 

9.2 Ensemble-Averaged Boundary Layer Velocity Distributions 

Boundary layer profiles were taken for one steady and two unsteady inlet flow 
conditions on the suction surface at 31 streamwise positions. Periodic unsteady flow was 
established by the wake generator, which included two clusters of cylindrical rods that were 
attached to the timing belts running with a translational speed of U = 5.0 m/s. To account 
for the unsteadiness caused by the frequency of the individual wake generating clusters and 
their spacing, the flow velocity, and the cascade parameters, the previously defined reduced 
frequency Q is used. The individual cluster configurations with the corresponding Q- 
parameter are specified in Table 4. 1 . Experimental investigations were perfonned for three 
different values of Q = 0.0, 1 .59, and 3.18. These values cover the entire reduced frequency 
range encountered in LPT-design and off-design operation conditions. 

Figure 9.3 displays two representative temporal ensemble-averaged velocity 
distributions for steady and two unsteady flow conditions with their characteristic features. 
All these figures show the boundary layer development from the core flow region to the 
blade surface at different streamwise positions. Approaching the wall surface, all velocities 
experience a continuous deceleration. The velocity gradient in both cases causes generation 
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and formation of vortices that transform steady nature of steady case into an unsteady one 
as clearly demonstrated in Figure 9.3. The unsteady case displayed in Figure 9.3 is 
characterized by its deterministic temporal periodicity. Approaching the wall surface, the 
traveling wake experiences a phase shift, while maintaining its deterministic nature. 
However, by penetrating into the boundary layer, the interaction between wake and 
boundary layer causes the detenninistic nature to degenerate into a stochastic one. As the 
flow continues downstream, continuous wake defonnation occurs in large velocity peaks 
inside the boundary layer. Convecting downstream, passing the separation zone, the wake 
diminishes, while the turbulent boundary layer starting from s/s o =0.766 dominates the flow 
picture. 
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Figure 9.3: Ensemble-averaged velocity as a function of time for three different 
reduced frequency of Q=0, 1.59, 3.18 at s/s o =0.0208, Re=l 10,000 
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Figure 9.5: Ensemble-averaged velocity as a function of time for three different 
reduced frequency of Q=0, 1.59, 3.18 at s/s o =0.546, Re=l 10,000 
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Figure 9.9: Ensemble-averaged velocity as a function of time for three different 
reduced frequency of Q=0, 1.59, 3.18 at s/s o =0.705, Re=l 10,000 
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Figure 9.10: Ensemble-averaged velocity as a function of time for three different 
reduced frequency of Q=0, 1.59, 3.18 at s/s o =0.766, Re=l 10,000 
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Figure 9.1 1: Ensemble-averaged velocity as a function of time for three different 
reduced frequency ofQ=0, 1.59, 3.18 at s/s o =0.951, Re=l 10,000 
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9.3 Time Averaged Velocity Distributions 

Consistent with the surface pressure distribution which is discussed above, the effect 
of the wake frequency on time-averaged velocity profiles and turbulence fluctuation 
distribution are presented for the corresponding rod spacings S R =°° mm, 80 mm, and 160 
mm at different longitudinal locations in Figures 9.12 to9.21. Upstream of the separation 
zone, the flow is attached to the blade, and the separation zone start at s/s o =0.52. Upstream 
of the separation zone at s/s o =0.52 and also at its immediate proximity s/s o =0.588, the 
velocity distributions inside the boundary layer experience a moderate decrease with 
increasing the reduced frequency. Inside the separation zone at s/s o =0.705, a substantial 
influence of the wake frequency is observed. Higher wake frequency introduces fluctuation 
kinetic energy into the boundary layer trying to reverse the separation tendency. And also, 
becalmed regions generated by the turbulent spots produced in the wake paths, were 
effective in suppressing or reducing the separation zone. As it can be seen from the velocity 
distribution profiles, the onset and the length of the separation zone is not changed. But, 
there is a slight change of the bubble height. This shows that flow does not have the 
capability to suppress the separation zone. It only reduces the separation bubble height. 
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Figure 9.12: Distribution of the time-averaged velocity along the suction surface 
for three different reduced frequencies of Q=0 (S R =°°), Q=1.59 (S R =160 mm), 
Q=3 . 1 8 (S r =80 mm) at Re=l 10,000 
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Figure 9.13: Distribution of the time-averaged velocity along the suction surface 
for three different reduced frequencies of Q=0 (S R =°°), Q=1.59 (S R =160 mm), 
Q=3 . 1 8 (S r =80 mm) at Re=l 10,000 
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Figure 9.14: Time-averaged velocity profiles along the suction surface of the blade at Q=0, Re=l 10,000 
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Figure9. 15: Time-averaged velocity profiles along the suction surface of the blade at Q=1.59, Re=l 10,000 
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Figure 9.16: Time-averaged velocity profiles along the suction surface of the blade at Q=3.18, Re=l 10,000 



Figure 9.17: Distribution of the time-averaged turbulence intensity along the 
suction surface for three different reduced frequencies of Q=0 (S R =°°), Q=1.59 
(S r =160 mm), Q=3. 18 (S R =80 mm) at Re=l 10,000 
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Figure 9.18: Distribution of the time-averaged turbulence intensity along the 
suction surface for three different reduced frequencies of Q=0 (S R =°°), 
Q=1.59 (S r =160 mm), Q=3.18 (S R =80 mm) at Re=l 10,000 
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Figure 9.19: Time-averaged turbulence-intensity profiles along the suction surface of the blade at Q=0, 
Re= 110,000 
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Figure 9.20: Time-averaged turbulence-intensity profiles along the suction surface of the blade at Q=1.59, 
Re= 110,000 
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Figure 9.21: Time-averaged turbulence-intensity profiles along the suction surface of the blade at Q=3.18, 
Re= 110,000 


Downstream of the separation zone, where the flow is fully reattached, s/s o =0.951, the 
impact of wake on the boundary layer is reduced. This effect is clearly shown in the velocity 
distribution at s/s o =0.951 . In accord with the previous investigations by Schobeiri et al. 
( 1 997) on a HP -turbine cascade, an increased wake frequency causes turbulence fluctuations 
to rise inside and outside the boundary layer as shown in Fig 9.17 to 9.21. 

9.4 Temporal Behavior of the Separation Zone Under Unsteady Wake Flow 

Velocity distributions on the suction surface with time as parameter are plotted in 
Figure 9.22 and 9.23. The non-dimensional time (t/x) values are chosen so that they 
represent the temporal states within one full period of wake passing. As Figures 9.22 (a) to 
(h) shows, the velocity distributions inside and outside the boundary layer at fixed s/s 0 - 
locations experience moderate to pronounced changes at Q=1.59. Starting with the initial 
state of flow being under the influence of free-stream turbulence to the final state where the 
flow is totally under the influence of wake the boundary layer undergoes a sequence of 
changes. Figure 9.22 (a) represents the instantaneous velocity distribution upstream of the 
separation zone followed by Figure 9.22 (b,c,d,e) which represent the velocity distributions 
inside the separation zone. Figure 9.22 (f) exhibits the instantaneous velocity distribution 
downstream of the separation zone. Figure 9.22 (g,h) represent the instantaneous velocity 
distributions, where the flow is fully reattached. In discussing the following results, we 
simultaneously resort to the wake distribution as well as the turbulence fluctuation results. 

Figure 9.22 (a) exhibits the velocity distribution on the suction surface at s/s 0 = 
0.402. At this streamwise position, the laminar boundary layer is subjected to a strong 
negative pressure gradient. The boundary layer distributions at different (t/x) experience 
changes in magnitude that reflect the corresponding changes of the impinging periodic wake 
velocity. It is worth noting, that despite the injection of turbulence kinetic energy by the 
vortical core of the impinging wakes, no local instantaneous boundary layer transition 
occurs. This is because of the strong negative pressure gradient that prevents the boundary 
layer from becoming instantaneously transitional. 

Representative for the instantaneous velocity distribution inside the separation zone, 
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shown in Figure 9.22 (b,c,d,e,f), we discuss the results plotted in Figure 9.22 (e ) at s/s 0 = 
0.674. During the time interval fromt/T close to 0.5 (1.5, 2.5, etc) to about t/x = 0.75 (1.75, 
2.75 etc) the separation zone is exposed to the wake external flow being under the influence 
of relatively lower turbulence does not have the capability to suppress the separation zone. 
Thus, the separation region is clearly shown by the velocity distributions at t/x = 0.5 and t/x 
= 0.75 . As the wake passes over the blade at s/s 0 = 0.674 introducing high turbulence kinetic 
energy into the boundary layer, the boundary layer is energized causing the separation zone 
to partially reduce or disappear. This leads to an instantaneous re-attachment. This time 
interval corresponds to the case where the flow is completely under the influence of wake 
and correspondingly the re-attached velocity distribution assumes a turbulent profde 
characterized by the curves at t/x = 1.0, t/x = 0.05, and t/x = 0.25 shown in Figure 9.22(e). 
To underline this statement, the steady state velocity distribution at the same streamwise 
position is also plotted in Figure 9.22 (e) using full circles. It shows clearly the separated 
nature of the boundary layer which coincides with the instantaneous velocity profile at t/x 
= 0.5. Intermediate times reflect the gradual change between the separation and re- 
attachment as the flow is undergoing the influence of the oncoming wake. Moving to the 
trailing edge of the separation zone, at s/s 0 = 0.705, when the flow is completely under the 
influence of the wake corresponding to the curves at t/x = 1.0, t/x = 0.05, and t/x = 0.25 
shown in Figure 9.22 (f), a partial reduction in boundary layer thickness as a result of wake 
impingement is visible, however, the separation zone does not seem to disappear. But, the 
compression of the bubble height is very clear. 

Figure 9.22 (g,h) exhibits the velocity distributions on the suction surface 
respectively at s/s 0 = 0.849 and s/s 0 = 0.898 where the flow is fully reattached. At this 
streamwise position, the boundary layer distributions at different (t/x) experience changes 
in magnitude that reflect the corresponding changes of the impinging periodic wake 
velocity. Once, the boundary layer is reattached and the velocity distribution assumes a fully 
turbulent profile, no major changes are observed. 

Figures 9.23 (a) to (h) show, the velocity distributions inside and outside the 
boundary layer at fixed s/s 0 -locations experience moderate to pronounced changes at 
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Figure 9.22: Distribution of the ensemble-averaged velocity development along 
the suction surface for different s/s 0 with time t/x as parameter for Q=1.59 
(S r =160 mm), Re=l 10,000 
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Figure 9.23: Distribution of the ensemble-averaged velocity development along 
the suction surface for different s/s 0 with time t/x as parameter for Q=3. 18 
(S R = 80 mm), Re=l 10,000 
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Q=3. 18. Similar instantaneous velocity distribution is observed when operating at a reduced 
frequency of Q=3 .18. Compared with Q= 1 .59, a stronger suppression of the separation zone 
is noticed at Q=3 .18. When the wake passing frequency is increased over the blade, it is also 
increasing high turbulence kinetic energy. Therefore, boundary layer is energized 
comparatively more than the reduced frequency of Q=1 .59 and it causes the separation zone 
to partially reduce more. 

9.5 Temporal-Spatial Resolution of the Separation Zone 

The separation zone can be thought of as a curve that connects the velocity inflection 
points along the suction surface. Starting with a reduced frequency of Q = 1.59 (S R = 160 
mm) at s/s o =0.520 , Figure 9.24 (a) exhibits the start of the separation zone with a lateral 
extension of about 0.2 mm. The impingement of the periodic wake vortical core with the 
high turbulence intensity causes a local periodic contraction of the zone in lateral direction. 
Convecting downstream, the lateral extension increases, thus the contraction appears more 
pronounced, Figure 9.24 and 9.25 (b,c,d). Substantial contraction occurs toward the trailing 
edge of the separation zone as shown in Figure 9.26 (e,f). At this streamwise position, the 
separation zone starts to contract at t/x=1.25 and 2.25. This points coincides with the 
streamwise position of the velocity maximum, which exactly correspond to the position of 
the minimum fluctuation. Similar results are observed when operating at a reduced 
frequency of Q=3.18 which is shown in Figure 9.27 to 9.29 (a,b,c,d,e,f) . As we discussed 
earlier from the other results, increasing the reduced frequency help to reduce more 
separation boundary layer. As a result, it helps to decrease the losses coming from the large 
bubble. 

Increased passing frequency causes an increase in turbulence kinetic energy. This in 
turn energizes the boundary layer. This is seen by comparing the velocity distribution for 
Q=3.18 and Q=1.59. 
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Q=1.59,s/s o =0.52 



(b) Q=1 .59,s/s o =0.546 



Figure 9.24: Contour plot of ensemble-averaged velocity distribution showing the 
effect of periodic wakes on the separation zone at Q=1.59 (S R =160 mm) at 
(a) s/s o =0.52 and (b) s/s o =0.546, Re=l 10,000 
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Figure 9.25: Contour plot of ensemble-averaged velocity distribution showing the 
effect of periodic wakes on the separation zone at Q=1.59 (S R =160 mm) at 
(c) s/s o =0.588 and (d) s/s o =0.617, Re=l 10,000 
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Figure 9.26: Contour plot of ensemble-averaged velocity distribution showing the 
effect of periodic wakes on the separation zone at Q=1.59 (S R =1 60 mm) at 
(e) s/s o =0.651 and (f) s/s o =0.674, Re=l 10,000 
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Figure 9.27: Contour plot of ensemble-averaged velocity distribution showing the 
effect of periodic wakes on the separation zone at Q=3.18 (S R =80 mm) at 
(a) s/s o =0.52 and (b) s/s o =0.546, Re=l 10,000 
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Figure 9.28: Contour plot of ensemble-averaged velocity distribution showing the 
effect of periodic wakes on the separation zone at Q=3.18 (S R =80 mm) at 
(c) s/s o =0.588 and (d) s/s o =0.617, Re=l 10,000 
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Figure 9.29: Contour plot of ensemble-averaged velocity distribution showing the 
effect of periodic wakes on the separation zone at Q=3.18 (S R =80 mm) at 
(e) s/s o =0.651 and (f) s/s o =0.674, Re=l 10,000 
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9.6 Change of Separation Bubble Height Under Unsteady Wake Flow 

The shape of the bubble configuration can be identified from the time-averaged 
velocity distributions. The dividing streamline and the line of inflection points is found by 
taking second derivative and equating to zero. A similar bubble formation is found for three 
different frequencies as shown in Figure 9.30. The separation points, where the separation 
starts and ends, remain the same for the three different frequencies. The separation bubble 
height decreases with the increase of the reduced frequency. The separation zone is exposed 
to wake external flow being under the influence of relatively lower turbulence followed by 
a periodic sequence of wake impingements. When the wake passes over the blade 
introducing high turbulence kinetic energy into the boundary layer, the boundary layer is 
energized causing the separation zone to partially reduce. Thus, this flow does not have the 
capability to completely suppress the separation zone. The main parameters describing the 
separation bubble for the three test frequencies are presented in Table 9.1. 



s/s o 

Figure 9.30: Change of separation bubble height under the influence of 
different reduced frequencies of Q=0, 1.59, 3.18 (S R =°°, 160 mm, 80 mm), 
Re=l 10,000 
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Table 9.1: Parameters describing separation bubble for three different frequencies 


Parameters 

Q=0.0 

Q=1.59 

Q=3.18 

s s /s 0 

0.5 

0.5 

0.5 

S n/ S o 

0.705 

0.705 

0.705 

s/s 0 

0.746 

0.746 

0.746 

Umm) 

6.123 

5.818 

5.358 


9.7 Separation Bubble Behavior Under Wake Propagation 

The effect of periodic unsteady wakes on the onset and extent of the separation zone 
is shown in Figure 9.31 to 9.36 for two different frequencies, namely Q=1.59 and Q=3.18. 
The wake propagation for Q=1 .59 and Q=3 . 1 8 is analyzed and the value of t/T corresponds 
to the point in the cycle at which data acquisition system is triggered. During a rod passing 
period, the wake flow and the separation zone undergo a sequence of flow states which are 
not noticeably different when the unsteady data are time-averaged. Starting with a reduced 
frequency of Q=1 .59 at t/x=0.05, the wake hits the blade leading edge and there is no effect 
on the separation zone. As the wake convects downstream, at t/x=0.50, the wake passing 
over the blade introducing high turbulence kinetic energy into the boundary layer, the 
energized boundary layer is partially reduced or disappeared. Becalmed regions generated 
by the turbulent spot produced in the wake paths, were effective in suppressing or reducing 
the flow separation. At s/s o =0.705 the core region has slightly moved towards to the leading 
edge. At t/x=l , the full effect of the wake on the boundary layer can be seen before another 
wake appears and the bubble moves back to the original position. Once the wake starts to 
penetrate the separation zone, the turbulent spot produced in the wake paths caused total 
suppression at some streamwise positions. Similar results are observed when operating at 
a reduced frequency of Q=3. 18. The wake frequency exerts no influence on the position of 
the separation zone. 
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(a) 


Q=1.59,t/T=0.05 




Figure 9.31: Ensemble-averaged velocity contours along the suction surface 
for different s/s 0 with time (a) t/T=0.05 and (b) t/x=0.25 as parameter for 
Q= 1. 59 (S R = 160 mm), Re=l 10,000 
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Figure 9.32: Ensemble-averaged velocity contours along the suction surface 
for different s/s 0 with time (c) t/x=0.50 and (d) t/x=0.75 as parameter for 
Q=1.59 (S r =160 mm), Re=l 10,000 
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(e) Q=1 .59, t/x=1 



Figure 9.33: Ensemble-averaged velocity contours along the suction 
surface for different s/s 0 with time (e) t/T=l as parameter for Q=1.59 
(S r =160 mm), Re=l 10,000 


(a) 


Q=3.1 8, t/x=0.05 



s/s 0 


Figure 9.34: Ensemble-averaged velocity contours along the suction 
surface for different s/s 0 with time (a) t/x=0.05 as parameter for Q=3.18 (S R =80 
mm), Re=l 10,000 
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(b) £2=3.18, t/x=0.25 



Figure 9.35: Ensemble-averaged velocity contours along the suction surface 
for different s/s 0 with time (b) t/x=0.25 and (c) t/l=0.50 as parameter for 
Q=3.18 (S r =80 mm), Re=l 10,000 
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Q=3.18,t/T=0.75 



Q=3.18, t/x=1 



Figure 9.36: Ensemble-averaged velocity contours along the suction surface 
for different s/s 0 with time (d) t/l=0.75 and (e) t/T=l as parameter for 
Q=3.18 (S r =80 mm), Re=l 10,000 
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9.8 Contraction, Separation and Regeneration of the Separation Zone 

Studying the temporal distribution of the turbulence fluctuations along the 
longitudinal extension of the separation zone, it becomes clear that, in conjunction with the 
pressure gradient and periodic wakes, there is another crucial mechanism responsible for 
the lateral and longitudinal decrease of the separation zone. This mechanism constitutes a 
combination of the high turbulence fluctuation level and its gradient. It is the temporal 
gradient of the turbulence fluctuation, or more precisely the fluctuation acceleration 
d v rms /dt that provides higher momentum and energy transfer into the boundary layer 
energizing the separation zone and causing it to partially or entirely disappear. 

For better understanding this phenomenon, the wake and the fluctuation velocity is 
presented in Figure 9.37 for the streamwise position s/s o =0.65 1 . The fluctuation nns velocity 
is scaled up with the factor of 4 to make it clear and an arbitrary normal position of 
y=2.85mm is chosen which is sufficiently above the separation zone for a given normal 
distance. 

Figure 9.37 (a) shows two separate regions: (1) a wake vortical core, occupied by 
vortices that originate from the moving cylindrical rods and generate high turbulence 
fluctuations, and (2) a wake external region between the adjacent vortical cores with 
relatively low turbulence activities. The wake configuration is asymmetric as shown in 
Figure 9.38. As Figure 9.38 shows, the separation zone starts to contract at t/l=1.25. This 
point coincides with the streamwise position of the velocity maximum, which exactly 
corresponds to the position of the fluctuation minimum, as shown in Figure 9.37 (b). At this 
point, the fluctuation within the vortical core starts to increase while the velocity 
continuously decreases. This process continues until the end contraction at t/T=l .41 has 
been reached. Thereafter, the separation zone is subjected to a process of intensive exchange 
of momentum and energy that causes the separation to diminish, as shown in Figure 9.38. 
The process of separation contraction, suppression, and regeneration is summarized in 
Figure 9.37 (b). It shows more details of separation contraction and suppression. In this 
context, it is necessary to subdivide the vortical core shown in Figure 9.37 (b) into four 
distinct regions, separated by thick dashed lines, as presented in Figure 9.37 (b). Region (a) 
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(a) Q=1 ,59,s/s o =0.651 , y=2.85mm 



(b) Q=1 .59,s/s o =0.651 , y=2.85mm 



Figure 9.37: Separation, contraction, suppression and regeneration at Re=l 10,000 
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is characterized by the initial positive gradient of the fluctuation dv rms /df> 0 marked with 
an upward arrow. Region (b) represents the substantial part of the vortical core with an 
intense turbulence activity. Region (c) serves as a transition between region (b) and the 
relatively calm region (d) characterized by dv rms /dt< 0 . 


E 

E 

> 


Q=1 .59, s/s o =0.651 



Figure 9.38: Separation zone, definition of contraction begin, contraction end, 
suppression, and regeneration 


For an initial fluctuation gradient dv nns ldt> 0, the separation zone begins to contract at 
t/x=1.25. This initial gradient is crucial for initiating the contraction process. Once the 
maximum fluctuation velocity with the temporal gradient<9v, mi /dt = 0 att/x=1.41 is reached 
(region a), the process of energizing the separation zone, thereby preventing its regeneration 
(region b). Passing the transition region (c), the process of suppression continues until the 
end of the vortical region at t/x=2.0 is reached. At this point the external wake region with 
its low turbulence content arrives causing a regeneration of the separation zone, thus 
reversing the entire suppression process. 
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While the existence of higher turbulence fluctuations expressed in terms of higher 
turbulence intensity is well known for influencing the flow separation, its gradient is crucial 
in suppressing or preventing the onset and the extension of the separation zone. The 
fluctuation gradient is an inherent feature of the incoming periodic wake flow and does not 
exist in a statistically steady flow, that might have high turbulence intensity. 

The results clearly indicate that for the particular blade under investigation, one has 
to do with a massive separation on the suction surface. These observations in comparison 
with the steady state reference case suggest that, once a massive separation zone is 
manifested, its size can be significantly reduced by periodic wake impingement, but it 
cannot be completely removed. The results presented here are valid for blades with a similar 
pressure distribution discussed earlier. Since the onset and extent of the separation zone is 
uniquely associated with the pressure gradient, blades can be designed with less local 
adverse pressure gradient, whose onset can completely be suppressed by impinging wakes. 

9.9 Ensemble-Averaged Fluctuation RMS Velocity Distribution 

The temporal-spatial contours of the ensemble-averaged fluctuation rms velocity at 
three different lateral positions for two reduced frequencies are presented in Figures 9.39, 
9.40 and 9.41. As shown in Figure 9.39, the wakes periodically disturb the boundary layer 
with high turbulence intensity cores. As Figure 9.39 (a) shows, three wake passings are 
visible for Q=1.59 at a nonnal distance from the wall of y=1.065m. At s/s o =0.52, the 
visibility of the wake is vanished due to separation. As explained earlier the separation zone 
starts at s/s o =0.52 and extends up to s/s o =0.746, thus occupying more than 24% of the 
suction surface and forming a massive separation. It is perfectly matching the result obtained 
from the Figure 9.1 and Figure 9.2. At s/s o =0.746, the fluctuation rms velocity field in 
Figure 9.39 (a) shows the stagnant fluid region and the area of the rms maximum, which 
indicates the development of transition and reattachment. Increasing Q to 3.18 causes an 
earlier mixing of the impinging wakes that leads to a complete degeneration of the 
deterministic periodic flow into a stochastic turbulent flow. For the highest reduced 
frequency of Q=3.18, the flow is highly turbulent all along the blade surfaces. As it can be 
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seen from the Figure 9.40 (b), increasing the reduced frequency has not brought major 
changes, regarding the structure of the separation zone. The laminar boundary layer and the 
onset of the separation zone remain the same. 

As shown in Figure 9.40 at y = 10.1 mm and Q=1.59, the wakes are visible up to 
s/s o =0.676. Up to this particular location, the wake convects above the separation zone, 
however by moving downstream the same lateral position is within the separation zone. This 
explains the sharp separation line where the wakes are not visible anymore. This implies that 
the wakes here, is in the previous figure have not penetrated into the separation zone. As 
shown in Figure 9.41, farther from the wall at y=60.1 and Q=1.59, the wake passings are 
visible, because the wake convects outside the boundary layer and the separation zone. 

9.10 Boundary Layer Ensemble-Averaged Integral Quantities 

The integral parameters, such as momentum thickness and shape factor, are of 
particular interest to turbine designer, since they provide an accurate first estimation of the 
quality of the designed blade. The ensemble-averaged distributions of the momentum 
deficiency thickness and shape factor for the suction surface at different streamwise 
positions are shown in Figure 9.42 and 9.43 for Q = 1 .59 . The momentum thickness values 
are non-dimensionalized with respect to the value corresponding to the steady case with Q 
= 0. The period T represents the wake-passing period that is specific to the individual wake 
generating cluster, which is characterized by the Q - value under investigation. The results 
are shown for three wake passing periods. The periodic behavior of the ensemble-averaged 
momentum thickness over the entire suction surface as a result of the embedded periodic 
wake flow, is clearly visible from Figure 9.42. It represents the momentum thickness 
behavior at different s/s 0 -locations upstream, within and outside the separation zone. The 
relative momentum thickness distribution upstream of the separation zone in Figure 9.42 
(b) integrally exhibits a slight increase, whereas inside the zone shown in Figure 9.43 (c), 
a substantial decrease is apparent. Atthe immediate vicinity of the trailing edge, s/s 0 = 0.705 
close to re-attachment, Figure 9.43 (d), the momentum thickness experiences a noticeable 
increase, which by convecting downstream decreases again and approaches the integral 
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(a) £1=1.59 



Figure 9.39: Ensemble-averaged rms fluctuation velocity in the temporal-spatial 
domain at y=1.065 mm for (a) Q=1.59 (S R =160 mm), and (b) 3.18 
(S r =80 mm), Re=l 10,000 


NASA/CR— 2009-21483 1 


83 


(a) £1=1.59 



Figure 9.40: Ensemble-averaged rms fluctuation velocity in the temporal-spatial 
domain at y = 10.1 mm for (a) Q=1.59 (S R =160 mm), and (b) 3.18 (S R =80 mm), 
Re= 110,000 
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(b) £1=3.18 
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s/s 0 


Figure 9.41: Ensemble-averaged rms fluctuation velocity in the temporal-spatial 
domain at y=60.1 mm for (a) Q=1.59 (S R =160 mm), and (b) 3.18 (S R =80 mm), 
Re= 110,000 
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values that are close to the steady state case. 

The ensemble averaged relative shape factor H 12 -distributions on the suction surface 
at the same streamwise locations are plotted in Figure 9.44 and 9.45 (a,b,c,d). Upstream of 
the separation zone, Figure 9.44 (b) they experience a similar periodic change with an 
average slightly close the steady case. Moving into the separation zone, Figure 9.44 and 9.45 
(b,c), each streamwise location presents its own form parameter that is specific to the 
velocity distribution. Similar results is observed when operating at a reduced frequency of 
Q=3.18 which is shown in Figure 9.46 to 9.49. 

9.11 Boundary Layer Time-Averaged Integral Quantities 

The time-averaged distribution of the boundary layer thickness, displacement 
thickness, the momentum deficiency thickness and the shape factor for the suction surface 
are shown in Figure 9.50 and 9.51. The momentum thickness values for three different 
reduced frequencies remain the same along the blade and do not reveal any unusual behavior 
within the separation bubble between the start s s and reattachment s r . The initial growth of 
the boundary layer displacement thickness and the shape factor, is comparatively small up 
to the maximum displacement location. As shown in Figures 9.50 and 9.51, both the 
boundary layer displacement thickness and shape factor are decreasing in the separation 
zone and after the reattachment as the reduced wake frequency is increased. This results 
clearly show the impact of the unsteady wake flow on the boundary layer parameters and 
hence on the profile loss coefficient and efficiency. Unsteady wake flow cause a reduction 
of the losses due to suppressed or reduced separated boundary layer. All of this results show 
that the steady state data cannot be transfonned to the unsteady turbine design technology 
without modification. 
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Figure 9.42: (a) and (b):Ensemble-averaged relative momentum thickness 
distribution along the suction surface for different streamwise positions 
at Q=1.59 (S r =160 mm), Re=l 10,000 
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Figure 9.43: (c) and (d):Ensemble-averaged relative momentum thickness 
distribution along the suction surface for different streamwise positions 
at Q=1.59 (S r =160 mm),Re=l 10,000 
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Figure 9.44: (a) and (b):Ensemble-averaged relative form parameter 
distribution along the suction surface for different streamwise positions 
at Q=1.59 (S r =160 mm) ,Re=l 10,000 
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Figure 9.45: (c) and (d):Ensemble-averaged relative form parameter 
distribution along the suction surface for different streamwise positions 
at Q=1.59 (S r =160 mm) ,Re=l 10,000 
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Figure 9.46: (a) and (b):Ensemble-averaged relative momentum thickness 
distribution along the suction surface for different streamwise positions 
at Q=3.18 (S r =80 mm) ,Re= 110,000 
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Figure 9.47: (c) and (d):Ensemble-averaged relative momentum thickness 
distribution along the suction surface for different streamwise positions 
at Q=3 . 1 8 (S r =80 mm) ,Re= 1 1 0,000 
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Figure 9.48:: (a) and (b):Ensemble-averaged relative form parameter 
distribution along the suction surface for different streamwise positions 
at Q=3 . 1 8 (S r =80 mm) ,Re= 1 1 0,000 
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Figure 9.49: (c) and (d):Ensemble-averaged relative form parameter 
distribution along the suction surface for different streamwise positions 
at Q=3 . 1 8 (S r =80 mm) ,Re= 1 1 0,000 
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Figure 9.50: Time-averaged boundary layer thickness and displacement thickness for 
three different reduced frequency of Q=0, 1.59, 3.18 ( S R =°°, 160 mm, 80 mm), 

Re= 110,000 
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Figure 9.5 1 : Time-averaged boundary layer momentum thickness and form 
parameter for three different reduced frequency of Q=0, 1.59, 3.18 ( S R =°°, 160 mm, 
80 mm), Re=l 10,000 
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9.12 Application of Boundary Layer Transition Models 

In order to predict the location of the boundary layer transition, some emprical 
transition models are used in CFD. One of them, is proposed by Mayle (1991) and 
Malkiel&Mayle ( 1 996) was included in this investigation. This model predicts the transition 
location based on the conditions at the separation point. A plot of the Reynolds number 
Re smd =LL(s md -s s )/v, based on the conditions at separation and the distance (s md -s s ) versus the 
momentum thickness Reynolds number Re 2s = U.,S 2 /v is shown in Figure 9.52. 

Re smd =700 Re 2s ° 7 for short bubbles 

Re smd =1400 Re 2s 0 ' 7 for long bubbles 

Short bubbles have only a local displacement. Therefore, the pressure distribution is close 
to that predicted value over the surface without separation. On the other hand, long bubbles 
interact with the exterior flow so, the pressure distribution over the blade surface is mostly 
different from the predicted one. Since long bubbles produce large losses and large 
deviations in exit flow angles, they should be avoided. On the other hand short bubbles are 
an effective way to force the flow to turbulent and may be considered as a means to control 
performance. The data are in good agreement with the empirical fonnula. 



Figure 9.52: Reynolds number based on the length of separation as a 
function of the momentum thickness Reynolds number at separation point 
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10. CONCLUSIONS 


A detailed experimental study on the behavior of the separation zone on the suction 
surface of a highly loaded LPT -blade under periodic unsteady wake flow was presented. One 
steady and two different unsteady inlet wake flow conditions with the corresponding passing 
frequencies, wake velocity and turbulence intensities were investigated utilizing a new large- 
scale, high-subsonic research facility The results of the unsteady boundary layer 
measurements were presented in ensemble-averaged, and contour plot forms. Surface 
pressure measurements were performed at Re= 50,000, 75,000, 10,0000, and 125,000. At 
each Reynolds number, one steady and two periodic unsteady inlet flow measurements were 
performed. The surface pressure distribution showed no major changes with respect to the 
above Re-number changes. Noticeable changes occurred, while operating at unsteady flow 
conditions. Detailed unsteady boundary layer measurement identified the onset and 
extension of the separation zone as well as its behavior under unsteady wake flow. Passing 
the wake flow with its highly turbulent vortical core over the separation region, caused a 
periodic contraction and expansion of the separation zone. It was found that, in conjunction 
with the pressure gradient and periodic wakes, the temporal gradient of the turbulence 
fluctuation, or more precisely the fluctuation acceleration dv /dt provides higher 
momentum and energy transfer into the boundary layer energizing the separation zone and 
causing it to partially or entirely disappear. We found that for dv rms /dt> 0 , the separation 
zone starts to contract whereas for dv rm Jdt < 0 it gradually assumes the shape before the 
contraction. The existence of higher turbulence fluctuations expressed in terms of higher 
turbulence intensity is well known for influencing the flow separation, its gradient is of 
crucial importance in suppressing or preventing the onset and the extension of the separation 
zone. The fluctuation gradient is an inherent feature of the incoming periodic wake flow and 
does not exist in a statistically steady flow, that might have high turbulence intensity. 
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A.l Data Reduction Programs 


1) Data reduction program for obtaining time-averaged parameter for steady 

measurement 

* 

* ENSEMBLE- This program calculates the time averaged velocity, 

* turbulence intensity, and fluctuation rms velocity 

* for single wire probes and applies ANGLE CORRECTION 

* FILES: flip : data file containing the the information on 

* (input) the number of data points and the names of 

* document file, data file and output file. 

* .dat -data files (input) 

* .doc -document filles. (input) 

* . ..a. out - output files, (output) 

* CTA.CAL : velocity calibration file (input) 

* Signal connection: 

* Chi : Trigger Ch3: Differential pressure transducer 

* Ch2 : Thermocouple Ch4: Hot wire 

* 

program ENSEMBLE 

parameter (max=2000000) 

integer nsamp,sampfreq,i,nb,j 

integer theta, theta 1 

real tsens,fcorr,vcorr,vcorrl 

real hwvel(4),sum,v 

real gamma, teal 

real trig, temp, pres, hw,kin_visc 

real ctc(5),acal(3),t(2) 

real u(max),du(max),turl(max),Re_ex(max) 
real dt(max),dx(max) 
character* 115 fil 
character* 115 inputfile 

c open the file containing velocity calibration constants for hot wire 
open(unit=8, file- eta. cal') 
read(8,*)(ctc(i), i=5,l,-l) 
write(*,*) ctc(5),ctc(4),ctc(3),ctc(2),ctc( 1) 
close(8) 
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open the file containing angle calibration coefficients 
open(unit=9,file='angcal.cof) 
read(9,*)(acal(i), i=l,3) 
close(unit=9) 
c open the master file 

open(unit= 1 0,file='BUSPROAC .PAR’) 

read(10,*) tsens 

read(10,*) theta 

read(10,*) teal 

close(lO) 

open(unit=13,file='caltherm.cof) 
read(13,*)(t(i), i=l,2) 
write(*,*) t(2),t(l) 
close(unit=13) 
open(unit=7,file='master.fil') 
read(7,*)nfile 

vcorr=acal(l)+(acal(2)*theta)+(acal(3)*theta*theta) 
vcorr 1 =acal( 1 )+(acal(2)*theta 1 )+(acal(3)*theta 1 *theta 1 ) 
fcorr=vcorr/vcorrl 

write(*,*)'V_correction factor is,',fcorr 

nsamp=20000 

sampfreq=20000 

do j=l,nfile 

read(7,200)inputfile 
do i=l,max 
u(i)=0.0 
enddo 

open(unit= 1 1 ,file=inputfile, status- old') 
sum=0.0 

do 510 i=l,nsamp 

read(l l,*)temp,hw 
temp=t(2) * temp+t( 1 ) 

hw=hw*sqrt((tsens-tcal)/abs(tsens -temp)) 
hw=hw * (hw * (hw * (hw * c tc(5 )+c tc (4))+ 

$ ctc(3))+ctc(2))+ctc(l) 

hw=4iw*fcorr 
hw=hw/vcorr 

c Calculate the sum for the ensemble average of the U velocity 
u(i)=hw+u(i) 
sum=sum+u(i) 

510 continue 
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close( 11) 
v=sum/nsamp 
do i=l,nsamp 

u(i)=u(i) 

enddo 

c Calculate the ensemble RMS and turbulence calculation 
do i=l,max 
du(i)=0.0 
enddo 

open(unit= 1 1 ,filc=inputfile, status- old') 
do 600 i=l,nsamp 

read(l l,*)temp,hw 
temp=t(2)*temp+t( 1 ) 
hw=hw * sqrt((tsens-tcal)/ abs(tsens -temp)) 
hw=hw * (hw* (hw * (hw * ctc(5 )+ctc(4))+ 

$ ctc(3))+ctc(2))+ctc(l) 

hw=hw*fcorr 
hwvel=hw/vcorr 
du(i)=(hw-v)**2 
600 continue 
close(l 1) 

delt=20 . 0/ float(nsamp) 
d=0.2715 
do i=l,nsamp 

du(i)=sqrt(du(i)) 
tur 1 (i)=du(i) * 1 00 ,/u(i) 
dt(i)=(i)*delt* 1000 
enddo 

ip=mdex(mputfile,' . ') 
fil=inputfile(:(ip-l))//'.out' 
c Write to the output fde 

open(unit=l 2, file=fil, status-new') 
write( 12,74) 
do i=l,nsamp 

write) 1 2,75)dt(i),u(i),du(i),turl (i) 
enddo 
close(12) 
enddo 
c Formats 
5 fonnat(a26,f9.3) 

74 fonnat(' dt[ms] u[m/s] du[m/s] Tu[%]’) 
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75 fomiat(4(lx,fl6.10)) 
200format(al 15) 
stop 
end 


2) Data reduction program for obtaining ensemble-averaged parameter for 
unsteady measurement 


ENSEMBLE- This program calculates the ensemble averaged velocity, 
turbulence intensity, and fluctuation rms velocity 
for single wire probes and applies ANGLE CORRECTION 
FILES: flip : data file containing the the information on 
(input) the number of data points and the names of 
document file, data file and output file. 

.dat -data files (input) 

.doc -document filles. (input) 

. ..a. out - output files, (output) 

CTA.CAL : velocity calibration file (input) 


Signal connection: 


Chi : Trigger 


Ch3: Differential pressure transducer 


* Ch2 : Thermocouple Ch4: Hot wire 
* 


program ENSEMBLE 

parameter (max=2500000) 

integer sampfreq,nsamp,i,nb,j 

integer nsamperevol,nrev 

integer ntext(6) 

real rps,rpm,vcorr,vcorrl,fcorr 

real hwvel(4) 

real gamma, tsens, teal 

real trig, temp, pres, hw,hwfil 

real ctc(5),acal(3),t(2) 

real u(max) ,du(max), tur(max) ,ufil(max) , dufil(max) 
real dt(max),dx(max),turfil(max) 
character* 50 fil 
character* 50 inputfile 

c open the file containing velocity calibration constants for hot wire 
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open(unit=8,file='cta.car) 
read(8,*)(ctc(i), i=5,l,-l) 
write(*,*) ctc(5),ctc(4),ctc(3),ctc(2),ctc(l) 
close(8) 

c open the file containing angle calibration coefficients 
open(unit=9,file='angcal.cof) 
read(9,*)(acal(i), i=l ,3) 
close(unit=9) 
c open the master file 

open(unit= 1 0,file='BUSPROAC .PAR’) 

read(10,*) tsens 

read(10,*) theta 

read(10,*) teal 

close(lO) 

open(unit=13,file='caltherm.cof) 
read(13,*)(t(i), i=l,2) 
write(*,*) t(2),t(l) 
close(unit=13) 
open(unit=7 , file- master . fil') 
read(7,*)nfile 

vcorr=acal(l)+(acal(2)*theta)+(acal(3)*theta*theta) 

vcorr 1=1. 6372625 

fcorr=vcorr/vcorrl 

write(*,*)'V_correction factor is = ',fcorr 
c write(*,*)acal(3),acal(2), acal(l) 
sampfreq=20000 
rps=l 
nrev=100 

nsamperevol=float(sampfreq)/rps 
nsamp=int(float(nrev)*nsamperevol) 
rpm=rps*60 
do j=l,nfile 

read(7,200)inputfile 
do i=l,max 
ufil(i)=0.0 
u(i)=0.0 
enddo 

open(unit= 1 1 ,file=inputfile,status='old') 
do 800 nb = l,nrev 

do 510 i=l,nsamperevol 
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c read(l l,*)temp,pres,hwfil,hw 

read( 1 1 ,*)trig,temp,pres,hwfd,hw 
temp=t(2)*temp+t( 1) 

hwfilHiwfil*sqrt((tsens-tcal)/(tsens -temp)) 
hwfil=hwfil*(hwfil*(hwfil*(hwfil*ctc(5)+ctc(4))+ 
$ ctc(3))+ctc(2))+ctc( 1 ) 

hw=hw*sqrt((tsens-tcal)/(tsens -temp)) 
hwHiw*(hw*(hw*(hw*ctc(5)+ctc(4))+ 

$ ctc(3))+ctc(2))+ctc( 1 ) 

hw=hw/fcorr 
c hw=hw/vcorr 

* Calculate the sum for the ensemble average of the U velocity 
ufil(i)=ufil(i)+hwfil 
u(i)=u(i)+hw 
510 continue 

800 continue 
close(ll) 

do i=l,nsamperevol 
ufil(i)=ufil(i)/nrev 
u(i)=u(i)/nrev 
enddo 

c Calculate the ensemble RMS and turbulence calculation 
do i=l,max 
du(i)=0.0 
dufd(i)=0.0 
enddo 

open(unit= 1 1 ,fde=inputfile, status-old') 
do 900 nb=l,nrev 
do 600 i=l,nsamperevol 

read(l l,*)temp,pres,hwfd,hw 
read) 1 1 ,*)trig, temp, pres, hwfd,hw 
temp=t(2) * temp+t( 1 ) 

hwfil=hwfil*sqrt((tsens-tcal)/(tsens -temp)) 
hwfilMrwfil*(hwfil*(hwfil*(hwfil*ctc(5)+ctc(4))+ 

$ ctc(3))+ctc(2))+ctc(l) 

hw=hw*sqrt((tsens-tcal)/(tsens -temp)) 
hw=hw * (hw * (hw * (hw * c tc(5 )+c tc (4))+ 

$ ctc(3))+ctc(2))+ctc(l) 

hw=hw/fcorr 
c hw=hw/vcorr 

dufd(i)=(hwfd-ufd(i))**2+dufd(i) 
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du(i)=(hw-u(i))* *2+du(i) 

600 continue 
900 continue 
close(l 1) 

delt=l .0/float(sampfreq) 
vbelt = 5.0*rpm/60.0 
delx = vbelt* delt 
do i=l,nsamperevol 

dufil(i)=sqrt(dufil(i)/nrev) 
turfil(i)=dufil(i)*100./ufil(i) 
du(i)=sqrt(du(i)/nrev) 
tur(i)=du(i)* 100./u(i) 
dt(i)=(i- 1 )* delt* 1 000 
dx(i)=(i-l)*delx*1000 
enddo 

ip=index(inputfde,'.') 
fil=inputfile( : (ip- 1 ))//' . out' 
c Write to the output fde 

open(unit=l 2, fde=fd, status-new') 
write( 12,74) 
write) 12,76) 
do i=l,nsamperevol 

write) 12, 75)dt(i),dx(i),ufd(i),dufd(i),turfd(i),u(i),du(i), 

& tur(i) 

c write) 12, 75)dt(i),ufd(i),dufd(i),turfd(i),u(i),du(i),tur(i) 

enddo 
close) 12) 
enddo 
c Formats 
5 fonnat(a26,f9.6) 

74 fonnat(’VARIABLES="dt(msec)", "dx(mm)", "ufil(m/s)", "dufil(m/s)", 
& ”Tufd(%)", "u(m/s)", ”du(m/s)", "Tu(%)'") 

76 fonnat(’ZONE 1=20000, F=POINT’) 

75 fonnat(8(lx,fl5.10) 

200 format(a50) 

stop 

end 


3) Data reduction program for separating the data into three segments and each 
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corresponding to one wake passing frequency 


* 

* SEPDATA - This program reads the data and separates 

* data into three separate segments, each corresponding 

* to one wake passing frequency 

* Files - Reads in LIMIT.TXT which specifies the starting and ending 

* length for each wake passing frequency. 

* Variables 

* nsamp=number of samples 

* sl=starting point (steady) 

* s2=starting point for the wake passing frequency zone (160 mm) 

* s3=starting point for the wake passing frequency zone (80 mm) 

* el=end point for the steady zone 

* e2=end point for the wake frequency zone (160 mm) 

* e3=end point for the wake frequency zone (80 mm) 


program sepdata 
implicit none 
integer max,ncount,nst,ip 
parameter(max=60000) 
integer nsamp,i,sampfreq,nfile,j 
real sl,s2,s3,s4,el,e2,e3,e4 
real dt(max) ,dx(max) ,u(max), du(max) , tur(max) 
real ufil(max), dufil(max) ,turfil(max) 
real dtm,delt 
character* 135 label 
character*60 fil 1 ,fil2,fil3 
character* 60 inputfile 
open(unit= 1, file-limit. txf, status- old') 
read(l,*)nsamp 
read(l,*)sampfreq !Hz 
read(l,*)sl,el 
read(l,*)s2,e2 
read(l,*)s3,e3 
c read(l,*)s4,e4 
close(l) 

open(unit=7, file-master 1 .fil') 
read(7,*)nfile 
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do j=l,nfile 

read(7,'(a)')inputfile 
delt=l ./float(sampfreq) 
ncount = 1 

ip=index(inputfile,'.') 

fil 1 =inputfile( : (ip- 1 ))//'a. out' 

fil2=inputfile( : (ip- 1 ))//'b . out' 

fil3=inputfile( : (ip- 1 ))//'c . out' 

open(unit=2,file=inputfile, status-old') 

open(unit=3 ,file=fil 1 ) 

open(unit=4,file=fil2) 

open(unit=5,file=fil3) 

read(2,'(a)')label 

read(2,'(a)')label 

do i=l,nsamp 

read(2,24)dt(i),dx(i),ufil(i),dufil(i),turfil(i),u(i),du(i),tur(i) 
c selecting the first wake passing frequency zone (steady) 
if((dx(i).ge.sl).and.(dx(i).le.el)) then 
if(ncount.eq. 1) then 
nst = i 
ncount = 2 
endif 

dtm = (i-nst)*delt* 1000+75 

write(3,24)dtm,dx(i),ufil(i),dufil(i),turfil(i),u(i),du(i),tur(i) 

endif 

c selecting the second wake passing frequency zone (160 mm) 
if((dx(i).ge.s2).and.(dx(i).le.e2)) then 
if(ncount.eq.2) then 
nst = i 
ncount = 3 
endif 

dtm = (i-nst)*delt* 1000+400 

write(4,24)dtm,dx(i),ufil(i),dufil(i),turfil(i),u(i),du(i),tur(i) 

endif 

c selecting the third wake passing frequency zone (80 mm) 
if((dx(i).ge.s3).and.(dx(i).le.e3)) then 
if(ncount.eq.3) then 
nst = i 
ncount = 4 
endif 

dtm = (i-nst)*delt* 1000+700 
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c dtm = (i-nst)*delt* 1000+700 

write(5,24)dtm,dx(i),ufil(i),dufd(i),turfil(i),u(i),du(i),tur(i) 
endif 
end do 
close(2) 
close(3) 
close(4) 
close(5) 
enddo 

24 fonnat(8(lx,fl5.10)) 
end 

4) Data reduction program for obtaining time-averaged quantities 


program average 

implicit none 

character* 85 outfile 

character* 65 inputfile 

real dt, dx , u,du, Tu, melt, m_dx , m u 

real m_du,m_Tu 

integer i,j,k,nfde, counter 

open( 1 , file- media2. inf ',status='old') 

read(l,*)nfde 

read( 1 ,'(a)')outfile 

opcn( 3, filc=outfilc, status-unknown') 
do i=l,nfile 
m_u=0.0 
m_du=0.0 
m_Tu=0.0 
read( 1 ,'(a)')inputfde 
opcn(2,nic=inputnie,status='old') 
read(2,*) 
counter=0 

do while (.not.eof(2)) 
read(2,*)dt,u,du,Tu 
counter=counter+ 1 
m_u=u+m_u 
m_du=du+m_du 
m_T u=T u+m_T u 
end do 
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m_u=m_u/ float(counter) 
m_du=m_du/ float(c ounter) 
m_T u=m_T u / flo at(c ounter) 
write(3 , 1 00)m_u,m_du,m_T u 
close(2) 
enddo 

1 00format(3 ( 1 x,f 1 8 . 1 0)) 
end 

5) Data reduction program for obtaining time-averaged boundary layer parameters 


c This program generates the ensemble averaged boundary layer parameters, 
c delta, delta l,delta2,delta3,H12,H23, etc. as a function of time (t/T : 0-3) 


program ENSBL2 

implicit none 

integer i,j,n,k, counter 

integer npol, start, ns, ks 

character*65 inputfile,fol 

real y,srod,ycoord,itime 

real dt,dx,ufil,dufil,tufil,t,tau,ul ,du,tu 

real a,g,u,fl,f2,f3 

real yintpol,uintpol,upot,ypot 

real break0,cscoef0, break l,cscoefl 

real break2,cscoef2,break3,cscoef3 

real low, high, ue, delta, delta 1 ,delta2,delta3 

real H12,H23,xt 

real* 8 acoeff(20),var(20),coeff(20) 

real csitg,csval 

real Udelta,S 1 ,y 1 ,y2,Udelta 1 

dimension dt( 1 00),ufil( 1 00),dufil( 1 00) 

dimension tufil(lOO) 

dimension u 1 ( 1 00),du( 1 00),tu( 1 00) 

dimension ypot(50),upot(50) 

dimension a( 1 00), g( 1 00), y( 1 00), u( 1 00), f 1 ( 1 00),f2( 1 00),f3( 1 00) 
dimension yintpol(500),uintpol(500) 

dimension break0( 1 00),cscoef0(4, 1 00),breakl ( 1 00),cscoef 1 (4,100) 
dimension break2( 1 00),cscoef2(4, 1 00),break3 ( 1 00),cscoef3(4, 1 00) 
c csakm Gives spline coefficients of the data points 
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c csitg Uses the spline coefficients to calculate the integral 
c iwkin Assigns larger workspace for the program 
c csval Gives the interpolated values using the spline coeffs 
c lsfit Fits the potential flow with a first order polynomial 

external csakm, csitg, iwkin, csval, lsfit 
call iwkin(43338) 

write(*,*)'Enter the approximate BL thickness in mm. ' 
read(*,*)ycoord 
ycoord=ycoord/ 1 000. 
open(unit=2, file- rnaster.fil', status-old') 
read(2,25) fol 
read(2,25) inputfile 
25 fonnat(a65) 
close(unit=2) 
c Input of the data 

open(unit=8,file=inputfile, status-old') 

read(8,*) 

read(8,*) 

i=0 

do while(.not.eof(8)) 
i=i+l 

read(8,*) a(i),ufil(i),dufil(i),tufil(i),ul(i),du(i),tu(i) 
a(i)=a(i)/l 000.0 
g(i)=ul(i) 
enddo 

closc(unit=8) 

open(unit=14,file=fo 1 ,status='new') 
write( 14,74) 

c Calculation of the potential flow fit 

c Initializing the parameters 

n=i-l 
start=0 
i=l 

do while(a(i).le.ycoord) 
if(a(i).lt.ycoord) then 
start=i+ 1 

elseif(a(i).eq.ycoord) then 

start=i 

endif 
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i=i+l 

enddo 

c Assigning the data to fit for potential flow 
do i=start,n 

ypot(i-start+ 1 )=a(i) 
upot(i-start+ 1 )=g(i) 
enddo 

c Fitting the potential flow 
npol=l 
ns=n-start+l 

call LSFIT(ns,ypot,upot,npol,acoeff,var,coeff,ks) 
c Calculating Ue and u(i) where u-inviscid=a(i)*coeff(2)+coeff(l) 
c fl(i) function to be integrated for displacement thickness 

c f2(i) function to be integrated for momentum thickness 

c f3(i) function to be integrated for energy thickness 

Ue=coeff(l) 
do i=l,n 
y(i)=a(i) 

u(i)=g(i) - (a(i)*coeff(2)+coeff(l)) + Ue 
fl(i)=l-u(i)/Ue 
f2(i)=u(i)/Ue*(l .0-u(i)/Ue) 
f3(i)=u(i)/Ue*(l .0-u(i)*u(i)/Ue/Ue) 
enddo 

c Calculation of the boundary layer thickness, delta using the 
c calculated velocity "u". 

c The boundary layer thickness is located using interpolated values 
c of the calculated velocity "u". 
xt=0.0 

call csakm(n,y,u,breakO,cscoefO) 
do i=l ,500 

xt=xt+y(n)/500 

yintpol(i)=xt 

uintpol(i)=csval(xt,n- 1 ,breakO,cscoefO) 
enddo 
i=l 

do while((abs(Ue-uintpol(i)).ge.0.01*Ue).and.(i.le.500)) 
i=i+l 
enddo 
k=i-l 

delta=yintpol(k) 
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c Limits of integration 
low=0.0 
high=delta 

call csakm(n,y ,f 1 ,break 1 ,cscoef 1 ) 
delta 1 =csitg(low,high,n- 1 ,break 1 ,cscoef 1 ) 
call csakm(n,y,f2,break2,cscoef2) 
delta2=csitg(low,high,n- 1 ,break2,cscoef2) 
call csakm(n,y,f3,break3,cscoef3) 
delta3=csitg(low,high,n- 1 ,break3 ,cscoef3) 

H 1 2=delta 1 /delta2 
H23=delta2/delta3 
delta=delta*1000. 
delta 1 =delta 1*1000. 
delta2=delta2* 1000. 
delta3=delta3* 1000. 

c S 1 is the slope of the potential core of u vs y 
c write(*,*)coeff(l),coeff(2),coeff(3),Ue 
Sl=coeff(2) 

c Find delta and Udel at y/delta=l .0 
c Up(i)=Sl*a(i)+Ue 
do i=l,n 

yl=a(i)* 1000.0 

y2=a(i+l)* 1000.0 if((yl.lt.delta).and.(y2.gt.delta))then 

c Udeltal=Sl*delta/1000+Ue 

Udelta=(delta-y 1 ) * (g(i+ 1 )-g(i))/i (y2-y 1 )+g(i) 

endif 

enddo 

c write(*,*)Udelta,Udeltal 

c Output data to the file "bounlay.par" to be read by the program law of the wall 
write( 1 4,34)delta/ 1 000, delta 1/1 000,delta2/ 1 000,delta3/ 1 000,Udelta, S 1 ,Ue 
34 fonnat(7(lx,fl5.10)) 

74 fonnatC delta delta 1 delta2 delt 

&a3 Udelta SI Ue’) 

stop 
end 

SUBROUTINE LSFIT(N,X,Y,NPOL, A, VAR, Coeff,KS) 

C NOTE: This file contains LSFIT, LSGET, and SIMQ 

C LSFIT IS A LEAST-SQUARES POLYNOMIAL FIT OF Y VS. X. THE 

POLYNOMIAL 
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C IS OF ORDER NPOL. 

C N=NUMBER OF DATA POINTS 
C X=INDEPENDENT VARIABLE 
C Y=DEPENDENT VARIABLE 
C NPOL=ORDER OF POLYNOMIAL TO BE FIT 

C A IS THE COEFFICIENT MATRIX WHICH MUST BE DIMENSIONED TO AT 
LEAST 

C (NPOL+ 1 ) * (NPOL+ 1 ) IN THE CALLING PROGRAM 

C VAR IS A DUMMY ARRAY WHICH MUST BE DIMENSIONED TO AT LEAST 
C NPOL+NPOL+ 1 IN THE CALLING PROGRAM. 

C Coeff= ANSWER VECTOR, I.E. THE COEFFICIENTS OF THE POLYNOMIAL. 
Coeff 

C MUST BE DIMENSIONED TO AT LEAST NPOL+NPOL+1 IN THE CALLING 
PROGRAM 

IMPLICIT REAL *8 (A-H,0-Z) 

REAL *8 A(l),Coeff(l),VAR(l) 

REALM X(1),Y(1) 

NPOLP 1 =NPOL+ 1 
NPOL2=NPOL+NPOL 
DO 60 I=0,NPOL2 
IP=I+1 
XAV=0.d0 
XYAV=0.d0 
DO 50 J=1,N 

XINC=l.dO 
if(i.ne.O) then 
DO 40 K=1,I 

XINC=XINC*X(J) 

40 CONTINUE 

end if 

X A V=X A V +XIN C 
X Y A V=X Y A V +XIN C * Y( J) 

50 CONTINUE 

V AR(IP)=X A V /N 
Coeff(IP)=XYAV/N 
60 CONTINUE 
ICOUNT=0 
DO 80 1=1, NPOLP 1 

DO 70 J=l, NPOLP 1 

ICOUNT=ICOUNT+ 1 
INDEX=I+ J- 1 
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A(IC OUNT)= V AR(INDEX) 

70 CONTINUE 

80 CONTINUE 

CALL SIMQ(A,Coeff,NPOLP 1 ,KS) 

RETURN 

END 

C 

SUBROUTINE LSGET(NPOL,Coeff,XINPUT,YOUT) 

C LSGET IS THE LEAST-SQUARES ’GET’ SUBROUTINE TO BE USED IN 
C CONJUNCTION WITH LSFIT. 

C NPOL = ORDER OF POLYNOMIAL 
C XINPUT = X POINT WHERE Y IS WANTED 
C YOUT = FITTED VALUE OF THE POLYNOMIAL AT X 
REAL *8 Coeff(l) 

REAL *4 XINPUT, YOUT 
NPOLP 1 =NPOL+ 1 
YOUT=0. 

XINC=1. 

DO 10 1=1, NPOLP 1 

IF (I.GT.l) XIN C=XIN C *XINPUT 
Y OUT=Y OUT +XINC * Coeff(I) 

10 CONTINUE 
RETURN 
END 

C 

C SUBROUTINE SIMQ 
C PURPOSE 

C OBTAIN SOLUTION OF A SET OF SIMULTANEOUS LINEAR EQUATIONS, 
C AX=B 
C USAGE 

C CALL SIMQ(A,B,N,KS) 

C DESCRIPTION OF PARAMETERS 

C A - MATRIX OF COEFFICIENTS STORED COLUMNWISE. THESE ARE 
C DESTROYED IN THE COMPUTATION. THE SIZE OF MATRIX A IS 
C N BY N. 

C B - VECTOR OF ORIGINAL CONSTANTS (LENGTH N). THESE ARE 
C REPLACED BY FINAL SOLUTION VALUES, VECTOR X. 

C N - NUMBER OF EQUATIONS AND VARIABLES. N MUST BE .GT. ONE. 
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C KS - OUTPUT DIGIT 
C 0 FOR A NORMAL SOLUTION 
C 1 FOR A SINGULAR SET OF EQUATIONS 
C REMARKS 

C MATRIX A MUST BE GENERAL. 

C IF MATRIX IS SINGULAR , SOLUTION VALUES ARE MEANINGLESS. 

C AN ALTERNATIVE SOLUTION MAY BE OBTAINED BY USING MATRIX 
C INVERSION (MINV) AND MATRIX PRODUCT (GMPRD). 

C SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED 
C NONE 
C METHOD 

C METHOD OF SOLUTION IS BY ELIMINATION USING LARGEST PIVOTAL 
C DIVISOR. EACH STAGE OF ELIMINATION CONSISTS OF INTERCHANGING 
C ROWS WHEN NECESSARY TO AVOID DIVISION BY ZERO OR SMALL 
C ELEMENTS. 

C THE FORWARD SOLUTION TO OBTAIN VARIABLE N IS DONE IN 
C N STAGES. THE BACK SOLUTION FOR THE OTHER VARIABLES IS 
C CALCULATED BY SUCCESSIVE SUBSTITUTIONS. FINAL SOLUTION 
C VALUES ARE DEVELOPED IN VECTOR B, WITH VARIABLE 1 IN B(l), 

C VARIABLE 2 IN B(2), , VARIABLE N IN B(N). 

C IF NO PIVOT CAN BE FOUND EXCEEDING A TOLERANCE OF 0.0, 

C THE MATRIX IS CONSIDERED SINGULAR AND KS IS SET TO 1 . THIS 
C TOLERANCE CAN BE MODIFIED BY REPLACING THE FIRST STATEMENT. 
C 

SUBROUTINE SIMQ(A,B,N,KS) 

IMPLICIT REAL *8 (A-H,0-Z) 

DIMENSION A(1),B(1) 

C FORWARD SOLUTION 
TOL=O.OdO 
KS=0 
JJ=-N 

DO 65 J=1,N 
JY=J+1 
JJ=JJ+N+1 
BIGA=0 
IT=JJ-J 
DO 30 I=J,N 

C SEARCH FOR MAXIMUM COEFFICIENT IN COLUMN 
IJ=IT+I 

IF ( AB S(BIGA)- AB S( A(I J))) 20,30,30 
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20 


BIGA=A(IJ) 

IMAX=I 
30 CONTINUE 

C TEST FOR PIVOT LESS THAN TOLERANCE (SINGULAR MATRIX) 
IF(ABS(BIGA)-TOL) 35,35,40 
35 KS=1 
RETURN 

C INTERCHANGE ROWS IF NECESSARY 
40 Il=J+N*(J-2) 

IT=IMAX-J 
DO 50 K=J,N 
I1=I1+N 
12=1 1+IT 
SAVE=A(I1) 

A(I1)=A(I2) 

A(I2)=SAVE 

C DIVIDE EQUATION BY LEADING COEFFICIENT 

5 0 A(1 1 )=A(1 1 )/BIGA 
SAVE=B(IMAX) 

B(IMAX)=B(J) 

B(J)=SAVE/BIGA 

C ELIMINATE NEXT VARIABLE 
IF(J-N) 55,70,55 
55IQS=N*(J-1) 

DO 65 IX=JY,N 

IXJ=IQS+IX 

IT=J-IX 

DO 60 JX=JY,N 
IX JX=N * ( JX- 1 )+IX 
JJX=IXJX+IT 

60A(IXJX)=A(IXJX)-(A(IXJ)*A(JJX)) 

6 5 B (IX)=B (IX) - (B ( J) * A (IX J)) 

C BACK SOLUTION 

70NY=N- 1 
IT=N*N 
DO 80 J=1,NY 
IA=IT-J 
IB=N-J 
IC=N 

DO 80 K=1,J 

B(IB)=B(IB)-A(IA) *B(IC) 
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IA=IA-N 
80 IC=IC-1 
RETURN 
END 


6) Data reduction program for obtaining equally distributed profile along the y axis 
by using Akima interpolation 


Program int 
implicit none 

integer i,nfils, ip ,j, counter ,k,nintv,n 
real srod,dt,dx, coord, uintpol 

real ufdintpol,duintpol,tuintpol,dufilintpol,tufdintpol 

real ufil,dufil,tufil,ul,du,tu 

real xt,newcoord 

real breakO,cscoefO,xdat,fdat 

real csval 

character*45 inputfile 
character* 85 outputfile 
dimension inputfde(50),outputfde(50) 
dimension dx(50,20000) 

dimension dt(50,20000),ufd(50, 20000), dufd(50, 20000) 
dimension tufd(50, 20000), coord(50),xt(50) 
dimension ul(50, 20000), du(50, 20000), tu(50, 20000) 
dimension ufdintpol(50, 20000), newcoord(50),xdat(50),fdat(50) 
dimension dufdintpol(50, 20000), uintpol(50, 20000), duintpol(50, 20000), 
& tuintpol(50, 20000), tufdintpol(50, 20000) 

dimension break0(50),cscoef0(4,50) 
c csakm Gives spline coefficients of the data points 
c iwkin Assigns larger workspace for the program 
c csval Gives the interpolated values using the spline coefficients 
external csakm, csval 
call iwkin(43338) 

opcn(unit=2, file-master, fif, status- old') 

read(2,*)nfds 

do i=l,nflls 

read(2,25) inpuffile(i) 

25 format(a45) 
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enddo 

close(2) 

open(unit=10,file- 1 l_234coord.dat', status-old') 
do i=l,nfils 

read(10,*) coord(i) 
enddo 
close(lO) 
do i=l,nfils 

open(unit=8,file=inputfde(i),status='old') 

j=0 

do while(.not.eof(8)) 
j=j+l 

read(8,*) dt(i,j),dx(i,j),ufil(i,j),dufil(i,j), 

& tufll(i,j),ul(i,j),du(i,j),tu(i,j) 

enddo 
counter=j 
enddo 

close(unit=8) 

n=nfils 

open(unit=20, file-10_56coord.dat', status-old') 
do i=l,nfils 

read(20,*) newcoord(i) 
enddo 
close(20) 
nintv=n- 1 

c Interpolation of the filtered velocity 
do j=l, counter 
do i=l,nfils 
fdat(i)=ufil(i,j) 
xdat(i)=coord(i) 
enddo 

call csakm(n,xdat,fdat,breakO,cscoefD) 
do i=l,nfils 

xt=newcoord(i) 

ufilintpol(i,j)=csval(xt,nintv,breakO,cscoefD) 

enddo 

enddo 

c Interpolation of the filtered nns velocity 
do j=l, counter 
do i=l,nfils 

fdat(i)=dufil(i,j) 
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xdat(i)=coord(i) 

enddo 

call csakm(n,xdat,fdat,breakO,cscoefO) 
do i=l,nfils 

xt=newcoord(i) 

dufdintpol(i,j)=csval(xt,nintv,breakO,cscoefO) 

enddo 

enddo 

c Interpolation of the filtered turbulence intensity 
do j=l, counter 
do i=l,nfils 

fdat(i)=tufil(i,j) 

xdat(i)=coord(i) 

enddo 

call csakm(n,xdat,fdat,breakO,cscoefD) 
do i=l,nfils 

xt=newcoord(i) 

tufilintpol(i,j)=csval(xt,nintv,breakO,cscoefD) 

enddo 

enddo 

c Interpolation of the raw velocity 
do j=l, counter 
do i=l,nfils 
fdat(i)=ul(i,j) 
xdat(i)=coord(i) 
enddo 

call csakm(n,xdat,fdat,breakO,cscoefD) 
do i=l,nfils 

xt=newcoord(i) 

uintpol(i,j)=csval(xt,nintv,breakO,cscoefD) 

enddo 

enddo 

c Interpolation of the raw rms velocity 
do j=l, counter 
do i=l,nfils 
fdat(i)=du(i,j) 
xdat(i)=coord(i) 
enddo 

call csakm(n,xdat,fdat,breakO,cscoefO) 
do i=l,nfils 
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xt=newcoord(i) 

duintpol(i,j)=csval(xt,nintv,breakO,cscoefO) 

enddo 

enddo 

c Interpolation of the raw velocity 
do j=l, counter 
do i=l,nfils 
fdat(i)=tu(i,j) 
xdat(i)=coord(i) 
enddo 

call csakm(n,xdat,fdat,breakO,cscoefD) 
do i=l,nfils 

xt=newcoord(i) 

tuintpol(i,j)=csval(xt,nintv,breakO,cscoefO) 

enddo 

enddo 

open(unit=22,file='masterout.fir, status-old') 
do i=l,nfils 

read(22,45) outputfile(i) 

45 format(a85) 
enddo 
close(22) 
do i=l,nfils 

open(unit=30,f"ile=outputfile(i), status-new') 
do j=l, counter 

write(30,24) dt(i,j),dx(i,j),ufilintpol(i,j),dufilintpol(i,j), 

& tufilintpol(i,j),uintpol(i,j),duintpol(i,j),tuintpol(i,j) 

24 fonnat(8( 1 x,f 1 5 . 1 0)) 

enddo 
enddo 
end 


7) Data reduction program for obtaining equally distributed values along y axis for 
time-averaged values 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
cThis program generates the interpolated time-averaged values along the boundary layer 
c 
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ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
program interpolation 
integer ndata,npoint 
parameter(ndata=5 0) 
parameter(npoint=28) 
integer i,nintv 

real break(ndata),cscoef(4,ndata),csval,f,xnew(ndata),fdata(npoint),xdata(npoint) 
real uintpol(ndata),tuintpol(ndata) 
real y(npoint),ufil(npoint),tufil(npoint),ul(npoint) 
real du(npoint),tu(npoint),dufil(npoint) 
character*40 inputfile 
character*30 outputfile 
external csakm,csval 
open(unit=2,file='master.fir,status='old') 
read(2,25) inputfile 
read(2,35) outputfile 
25 format(a40) 

35 fonnat(a30) 
close(unit=2) 

open(unit=8,file=inputfile,status='old') 
do i=l,npoint 

read(8,*)y(i),ul(i),tu(i) 

enddo 

close(unit=8) 

open(unit=9, file-newpoints.txt', status-old') 
do i=l,ndata 

read(9,*) xnew(i) 
enddo 

c Interpolation of the velocity 
do 10 i=l,npoint 
xdata(i)=y(i) 
fdata(i)=ul(i) 

10 continue 

call csakm(ndata, xdata,fdata, break, cscoef) 
nintv=ndata-l 
do 20 i=l, ndata 
x=xnew(i) 

uintpol(i)=csval(x,nintv, break, cscoef) 

20 continue 

c Interpolation of the turbulence intensity 
do 30 i=l,npoint 
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xdata(i)=y(i) 

fdata(i)=tu(i) 

30 continue 

call csakm(ndata, xdata,fdata, break, cscoef) 
nintv=ndata-l 
do 40 i=l, ndata 
x=xnew(i) 

tuintpol(i)=csval(x,nintv, break, cscoef) 

40 continue 

open(unit= 1 0 , fde=outputfde, status-new') 
do i=l, ndata 

write( 10,75) xnew(i),uintpol(i),tuintpol(i) 
enddo 

75 fonnat(3(lx,fl5.10)) 
end 

8) Data reduction program for obtaining contour plots from ensemble-averaged 
results 

c 

c Used for data analysis for unsteady boundary layer (UBL-series) 
c development study. 

c This is a processing program that converts s series of ensemble 
c averaged xxxxxxxa,b.OUT files to plots in the distance/time 
c domain. 


program CONT 

implicit none 

integer nfils,i,j,n 

real ufil,dufil,Tufil 

real dt,dx,u(50),du,Tu(50),y(50) 

real ratio 

character*30 inputfile 
character*35 outfile 
nfils=30 

outfile='D:\results\cfd\contourcfd.out' 
opcn(3,nie=outfile, status- unknown') 
write(3,74) 
write(3,76) 
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do i=l,50 

open(unit=l l,file='ratio.inf,status='old') 
do j=l,nfils 

open(unit=10,fde='contour.fd',status='old') 
read(l l,*)ratio 
read( 1 0,’(a)')inputfde 
open(2,file=inputfile,status='old') 
do n=l,50 

read(2 , * )y (n) ,u(n) ,Tu(n) 
enddo 

write(3,93)ratio,y(i),u(i),Tu(i) 
enddo 
close(lO) 
close(l 1) 
enddo 
close(3) 

93 fonnat(4(lx,fl5.10)) 

74 format('VARIABLES="s/s_o","y(mm)"," U (m/s)", M T U (%)'") 

76 fonnat(’ZONE 1=1500, F=POINT’) 
stop 
end 

9) Data reduction program for obtaining contour plots from time-averaged results 


c 

c Used for data analysis for unsteady boundary layer (UBL-series) 
c development study. 

c This is a processing program that converts s series of ensemble 
c averaged xxxxxxxa,b.OUT files to plots in the distance/time 
c domain. 


program CONT 
implicit none 
integer nfils,i,j 

real ufil,dufil,Tufil,m_Tufil,m_dufil,m_ufil 
real dt, dx , u,du, Tu, melt, m_dx , m u 
real ratio, srod, period, wallnonnal 
character* 75 inputfile,outfile 
write(*,*)'Enter rod spacing in mm.' 
read(*,*) srod 
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period=srod/5 

open(unit=10,fde='contour.fd',status='old') 

read(10,*)nfils 

read( 1 0,’(a)’)outfde 

open(3,fde=outfile,status=Tinknown') 

write(3,74) 

write(3,76) 

open(unit=l 1 , file-ratio, inf , status-old') 
c open(unit=l 1 ,file- 1 l_7coord.dat' ,status='old') 
do i=l,nfils 
read(l l,*)ratio 
c read( 1 1 ,*)wallnormal 
read( 1 0,’(a)')inputfile 
open(2,file=inputfile,status='old') 
do while (.not.eof(2)) 

read(2,*)dt,dx,ufil,dufil,Tufil,u,du,Tu 
c if(dt.ge.400)then 

if(dt.ge.700)then 
c dt=dt-400.0 

dt=dt-700.0 
if(dt/period. gt.3 . 0)then 
go to 22 
endif 

write(3,93)ratio,dt/period,ufil,dufil,Tufil,u,du,Tu 
c write(3,93)dt/period,wallnonnal,ufil,dufil,Tufil,u,du,Tu 

endif 
enddo 

22 close(2) 
enddo 
close(3) 

93 format(8( 1 x,f 1 5 . 1 0)) 

74 fonnat('VARIABLES="s/so","t/tau","ufrl(m/s)", M dufil(m/s)", "Tufil 
&(%)", "u(m/s)", M du(m/s)","Tu(%)"’) 
c74format('VARIABLES="t/tau","y/l","ufil(m/s)", M dufil(m/s)", "Tufil 
c &(%)", "u(m/s)","du(m/s)","Tu(%)"’) 

76 format('ZONE 1=29791, F=POINT') 
stop 
end 
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APPENDIX-B 
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B.l Uncertainty Analysis for Single Hot-Wire Measurements 
Kline-McKlintock Uncertainty Analysis for Single Hot-Wire Probe Calibration 
I) Calibration 


Tail- := (22.7+ 273. 15)K coTair := 0.5K 


Free-stream air temperature 


„ , . 753.8V 

Pamb := | in Hg 

25.4 ) ~ 


coPamb := | — \n Hg Atmospheric pressure 
25.4 ) ~ 


R := 287 


joule 

kg-K 


Gas constant, uncertainty assumed to be zero 


Ea := 0.33 


coEa := 0.002 


Voltage of Anemometer 


D1 := 38.1mm 


coDl := 0.5mm 


Outer diameter of calibration nozzle 


D2 := 152. Shim 


coD2 := 0.5mm 


Inner diameter of calibration nozzle 


( 0.04A 

Ptrans := in Hg 

1 25.4 ) ~ 


CoPtrans := 0.6Pa 


Pressure read from pressure transducei 


1) Determine uncertainty for air density 


P := 


Pamb 
R- Tan- 


ks 

p = 1.18— 
3 
m 


Note: DrhoDTair=Partial derivative of rho with respect to Tair 


DpDTair:= 


-Pamb 

RTair 2 


DpDTair= -3.989x 10 


3 kg 


m 3 K 


DpDPamb := 


1 


RTair 


DpDPamb = 1.1 78x 10 5 — 

2 

m 
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cop : 


2 

: (DpDTaircoTair) + (DpDPambcoPamb) 


-,0.5 


- 3 kg 

cop = 2.143x10 — 

3 
m 


%p := 


^jioo 

. p ) 


%p = 0.182 


V := 


(2- Ptrans ) 


n0.5 


D1 




D2 4 ) 


V = 3.008 — 


DVDD2 := -2.8284271247461 


Ptrans 


0.5 


0.5 


1 - 


D1 




1.5 


D2 4 ) 


D1 


D2 ' 


DVDD2 =-0.152 Hz 


DVDD1 := 2.8284271247461 


Ptrans 


0.5 


0.5 


1 - 


D1 




1.5 


D2 4 ) 


D1 


D2 


DVDD1 =0.611 Hz 


DVD p := -0.707106781186 •- 


Ptrans 


0.5 


1.5 


D1 




0.5 


DVD p = -1.274 


kg s 


D2 4 ) 


DVDPtrans= — 


0.7071067811865475244 


0.5 

0 

f 4A°- 5 ! 

D1 1 

r 

v D2 'J _ 


Ptrans 


0.5 


2 

DVDPtrans= 0.283 — - 
kg 
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coV := 


.l 0 - 5 

(DVDDlcoDl) + (DVDD2coD2) + (DVEpcop) + (DVDPt rans- coPtrans ) J 


m 

coV = 0.17— 


%V:= 


100 
V J 


%V = 5.642 


2 ) 
Determine 
uncertaint 

y for 

velocity 


The percentage of uncertainty for other velocities 


V := 15- 


%V:= 0.284 


m 

V := 10 — 


%V:= 0.529 


V:= 5 - 


%V:= 2.053 


%Vn:= 


toV- 


Lv 


lopy 

Vj 


-i0.5 



%Vn= 5.768 


The uncertainty of the curve fit according to S. Yavuzkurt (1984) is 1.2%. Thus the new 
velocity uncertainty is given by; 
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The curve fit uncertainty is calculated as follows: 

1 . The measured velocity is curved fit against anemometry voltage 

2. The residuals of the curve fit are calculated 

3. The residuals are divided by measured velocity at which they occur 

4. The quantity in 3 is squared and summed for all velocities 

5. The sum of the squares is divided by the number of data points and the square root 
is taken 
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6 . 


The percentage uncertainty in the curve fit is called beta and the equation is given by 


Determine uncertainty of velocity due to error in anemometer voltage 




2 3 zA m 

Vnl:= \1.374268+ 4.3283137Ea+ 1.4779551Ea - 0.24928452Ea + 0.1867769Ea )- 

s 


DVnlDEa:= 4.3283137+ 2.9559102Ea - 0.74785356Ea 2 + 0.7471076Ea 3 


DVnlDEa= 5.249 


coVnl := DVnlDEacoEa-- 


Vnl = 2.957- 


coVnl = 0. 01- 


Total calibration uncertainty for velocity due to velocity reading, curve fit, and voltage error 


%Vtotalc := 


%Vn 2 + 


100 


ooVnl Y 
Vnl ; J 


-,0.5 


%Vtotalc = 5.779 


S. Yavuzkurt is below 


II) Determine uncertainty of the velocity 


E:= 0.33 

coE := 0.002 

Tsen := 250 

coTsen := 0.5 

T :=24 

coT := 1 

Teal := 22.7 

coTcal := 1 


En is the anemometer voltage that has been 


after data reduction 

Anemometer voltage 
Temperature of hot-wire sensor (C) 
Temperature of air during measurements (( 
Temperature of air during calibration (C) 
corrected for temperature 
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(Tsen - T) 
(Tsen - Teal) 


En= 0.329 


DEnDTcal:= 0.5 E- 


(Tsen - T) " 
(Tsen - Teal) 1 


DEnDTcal= 7.238x 10 


DEnDT:=-0.5- 


(Tsen - T)°' 5 (Tsen - Tcal)°' 5 _ 


DEnDT= -7.28x 10 


DEnDE:= 


(Tsen - T) 
(Tsen - ’leal)* 


DEnDE = 0.997 


2 2 

COEn := (DEnDTcal coTcal) + (DEnDT coT) + (DEnDE- coE) 


CoEn = 2.243 x 10 


%En := 100- 


%En = 0.682 


Determine uncertainty of effective velocity due to error in corrected voltage 


Vn:=( 1.374268+ 4.3283137En+ 1.4779551En 2 - 0.24928452En 3 + 0.1867769En 4 ) — 

s 


Vn = 2.952- 


DVnDEn:= 4.3283137+ 2.9559102En - 0.74785356En 2 + 0.7471076En 3 


DVnDEn= 5.247 


COVdata := DVnDEncoEn - 


coVdata = 0.0 12- 


Total velocity uncertainty for calibration and data reduction 


f co V data 2 

%Vtotal := 100 + %Vtotalc 

LI Vn j 


%Vtotal = 5.792 
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